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The Unseen Population of F to K-type Companions to Hot 
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ABSTRACT 

We present a method to select hot subdwarf stars with A to M-type companions 
using photometric selection criteria. We cover a wide range in wavelength by combin- 
ing GALEX ultraviolet data, optical photometry from the SDSS and the Carlsberg 
Meridian telescope, near-infrared data from 2MASS and UKIDSS. We construct two 
complimentary samples, one by matching GALEX, CMC and 2MASS, as well as a 
smaller, but deeper, sample using GALEX, SDSS and UKIDSS. In both large 
number of composite subdwarf plus main-sequence star candidates were found. We 
fit their spectral energy distributions with a composite model in order to estimate the 
subdwarf and companion star effective temperatures along with the distance to each 
system. The distribution of subdwarf effective temperature was found to primarily lie 
in the 20, 000 — 30, 000 K regime, but we also find cooler subdwarf candidates, making 
up ~ 5 — 10 per cent. The most prevalent companion spectral types were seen to be 
main-sequence stars between FO and KO, while subdwarfs with M-type companions 
appear much rarer. This is clear observational confirmation that a very efficient first 
stable Roche-lobe overflow channel appears to produce a large number of subdwarfs 
with F to K-type companions. Our samples thus support the importance of binary 
evolution for subdwarf formation. 

Key words: stars: fundamental parameters - subdwarfs - white dwarfs - ultraviolet: 
stars - infrared: stars. 



1 INTRODUCTION 



Subluminous blue stars were first discovered by 
iHumason fc Zwickvl (|l947l ) in a photometr i c sur vey of 
the North Galactic Pole region. I Green et al l (|1986D found 
many more hot subdwarfs in the Palomar-Green (PG) 
survey, to the extent that they were the dominant species 
among faint (B < 16.1) blue objects. In the PG survey they 
outnumber white dwarfs (WD) and are prevalent enough 
to account for the u ltraviolet upturn in early-type galaxies 
ijBrown et alJ Il997l ). Hot subdwarf stars are either core 
helium-burning stars at the end of the horiz ontal branch 
or h a ve evolved even beyond that stage |Heber et al.l 
1 19841 : iHeberl 1 19861 ). They have a relatively well defined 
mass around the canonical (theoret i cal) value of 0.46 M 
jSaffer et all 1 19941 : lHan et all 120031 : iPolitano et"aH 1200; 
and radii of a few tenths of a solar radius. Their very thin 
layers of hydrogen (M env < O.OIM©) are not able to support 
shell burning after helium-core exhaustion. Thus instead of 



following the asymptotic giant branch route, they evolve 
more or less directly into white dwarfs. Observationally two 
classes are defined, those with helium-poor spectra (sdBs) 
and those that are helium-rich (sdOs). 

Formation scenarios of subdwarfs invoke either fine- 
tuned single star evolution or rely on close-binary star inter- 
actions. In the late hot-flasher scenario, a low-mass star un- 
dergoes the He core-flash at the tip of the red-giant branch. 
However, if sufficient mass is lost on the red giant branch, the 
star can experience the He c ore-flash whilst descending th e 
white dwarf cooling track (Cast ellani &: Castellanil Il993h . 
Such a star would end up close to the He main sequence 
(MS), at the very ho t end of the extreme horizontal branch 
i D'Cruz et al.| [l996h Alternatively, the formation involves 
one or two phases of common-envelope evolution and/or sta- 
ble Roche-lobe overflow (RLOF) within a close binary sys- 
tem (|Mengel et al.|[l976l h Binary evolution could even take 
the route o f merging two helium white dwarfs followed b y 
He ignition (|Webbink|| 19841 : Ilbenlll990l : ISaio fc JeffenfeoOCh . 
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All formation scenarios require substantial mass loss before 
the start of core He-burning, however the specific physical 
mechanisms for this are still unclear. A detailed rev iew on 
this and the field as a whole is given bv lHeberl <|2009h . 

Since the first quantitative estimates of the contribution 
of dif ferent binary channels to th e population of subdwarf 
stars (|Tutukov fc Yungelson|[l990l ). it has been shown that a 
large fraction of subdwarfs do reside in binaries. In the PG 
sample of subdwarfs, a significant fraction show compos- 
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percent: kllardetal. 


1 19941). Radial veloc- 


iMaxted et al. 2001; 


Morales-Rueda et al.l 



120031 ) confirm the high fraction of binaries with ra- 
tios as high as two-thirds. High-resolution optical spec- 
tra from the ESQ Sup ernova la Progenitor Survey (SPY; 
iNapiwotzki et all 120011), led to bina r y star fractions o f 
30-40 per cent ilNapiwotzki et all |2004| ; iLisker et all [200II) . 
ICopperwheat et al. ( 2011 ) estimate that the binary fraction 
in the sdB population is somewhat higher at 46 - 56 per cent. 
This is only a lower limit since the radial velocity variations 
that lCopperwheat et al.l (|201lT ) search for would be difficult 
to detect in long period systems. 

Other searches h a ve used near-infrared photometry 
(e.g iTheill et al.lll995l ; lUlla fc Thei3ll99St IWilliams et all 
l200lf ) or photometric catalogue s such as the Two Mi- 
cron All Sky Survey (2MASS; ISkrutskie et~aH 120061) to 



find subdwarfs w i th companion s (e.g . IStark fc Waddl2003l; 
iGreen et al.ll2006l ; IVennes et al.ll201ll ). Call absorption can 
also be used to infer the pre sence of a cooler companion 
star jjefferv fc Pollaccol [l998). The majority of compan- 
ions fo und to date have either been M-type stars or white 
dwarfs l|Hebej |2009). However, some F, G and K-type com- 
panions to subdwarfs hav e been seen in studies s uch a s 



Aznar Cuadrado fc Jeffervl (l200ll) [ Reed fc Stienind (|2004|) 
LiskeretalJ (120051) IWade et alj (12006b, IStark fc Wadd 
ll2006h.lWade et al.l (fioobl l. jMoni Bidin fc Piottol (|2010h and 
iGeier et all (MUCHFUSS; l2011bl ). Depending on the study, 
and its corresponding selection effects, the companions to 
subdwarfs have been shown to be mostly main-sequence 
stars (e.g. lAznar Cuadrado fc Jeff erv 2001, ) and giant or 
subgiant compa nions in some case s (e.g. lAllard et al.lll994l 
and BD-7°5977; lHeber et al.ll2002T l. 

Many of the previous surveys have be en biased by selec - 
tion effects and inhomogeneous data sets. lHan et alj (|2003l ) 
argued that a large number of sdB stars may be missing from 
current samples. Early-type main-sequence stars of spectral 
type A and earlier would outshine a subdwarf at optical 
wavelengths. F to K-type companions on the other hand, 
have generally been avoided because the spectral analysis of 
the composite spectrum becomes difficult. Systems with ear- 
lier type companions are actually predicted, in some cases, 
to be far more common than the M-type companions that 
have primarily been found so far. In the lHan et all (|2003l ) 
study, subdwarfs with early type companions are produced 
in the very efficient first stable RLOF channel and are ex- 
pected to be in syste ms with subdwarfs as cool as 15, 000 K. 
IClausen et all (|2012h , however, do not find the same multi- 
tude of F-type companions. Identifying this predicted pop- 
ulation, and determining their relative contribution to the 
subdwarf population would offer important constraints on 
the prior binary evolution that led to their formation. In 
addition, the distribution of orbital periods and subdwarf 



temperatures of such a sample will provide direct constraints 
on key param eters that underpin s ubdwarf population syn- 
thesis models l|Clausen et al.ll2012T ). 

In this study, we take advantage of recent large-area ul- 
traviolet, optical and infrared photometric surveys to search 
for new composite systems comprised of subdwarfs plus 
main-sequence star companions of mid-M-type and earlier. 
Cuts in colour-colour space are employed to separate these 
objects from possible contaminants. We also develop a fit- 
ting technique to simultaneously determine the subdwarf 
and companion effective temperatures from the photometric 
magnitudes. This permits the recovery of composite systems 
with much earlier type companions than seen in previous 
studies. Furthermore, we are sensitive to a wide range of 
separations and binary periods in that we only limit our- 
selves to spatially unresolved systems. Finally, we discuss 
the distribution of objects in effective temperature and dis- 
tance to the system. 



2 SYNTHETIC MODELS 

To aid our search for subdwarfs with companions, we pro- 
duced a grid of synthetic sdB and main-sequence star spec- 
tra, which allowed us to produce synthetic colours of the 
composite systems. 

The sdB spectra were c alculated using the model at- 
mosphere code described by iHeber et al.l ((2000), covering 
Tea = 11,000 - 40, 000 K in steps of 1,000K. The cor- 
responding surface gravities were chosen to ensure that 
our temperatu re sequence tracks th e (extreme) horizontal- 
branch stars (|Dorman et al. I Il993l ). This translates into 
log g = 4.0 for T eff = 11, 000-13, 000 K objects, log g = 4.5 
for T cff = 14, 000 - 16, 000 K, log g = 5.0 covering T cft = 
17, 000 - 20, 000 K, log g = 5.5 for T cff = 21, 000 - 28, 000 K 
and log g = 6.0 for T cff = 29, 000-40, 000 K. Surface gravity 
does not significantly affect spectral slope, but does affect 
the width of line profiles, which is a negligible feature when 
fitting photometry as we do here. It also corresponds to a 
significant change in the size of the subdwarf and therefore 
the relative brightness of the subdwarf and the companion. 

A range of solar metalicity main-sequence star tem- 
plates of effective temperatur es from 4, 250 K to 25, 000 K in 
48 steps were taken from the lCastelli fc Kuruczl (|2003l ) AT- 
LAS9 model a tmosphere library. For models below 4, 250 K, 
iPickled if 19981 ) stellar spectral libr ary models are subs t ituted 



because of the problems with Castelli & Kurucz (200 



model colours in this region l|Bertone et al.ll2004l) . A fPicklei 
(1998) MOV star is used as a proxy for a 4, 000 K model. Sim- 
ilarly, M1V, M2V, M3V and M5V replace 3, 750 K, 3, 500 K, 
3, 250 K and 3, 000 K models, respectively. We restrict the 
models to unevolved main-sequence stars because, as we 
will see in Section l4~2l sub-giant and giant companions do 
not contribute significantly to our sample. The impact of 
this will be further discussed in Sections]?] and [8] Both the 
sdB and main-sequence star spectra co ver the wavelength 
range 1, 150 — 25, 000 A. To normalise the lCastelli fc Kuruczl 
(2003) main-sequence star models to a flux at 10 pc, we 
rescale the models to match the l uminosities from the (zero 
age main sequence) isochrones of iGirardi et al. I (|2000h . 

The two grids of spectra were folded through all rel- 
evant filter transmission curves to calculate absolute mag- 
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nitudes. The two components could therefore be added at 
a common distance. To separate composite subdwarf plus 
companion systems from single subdwarfs and single main- 
sequence stars in colour-colour space, a large wavelength 
range must be sampled. The combination of a very blue and 
a red colour allows for a significant contribution from both 
the subdwarf and companion components to be seen in a 
colour-colour diagram. We therefore chose to cross-match an 
ultraviolet survey with a series of optical and near-infrared 
surveys. 

3 CROSS-MATCHING 

3.1 Sample I: GALEX, CMC and 2MASS 

The Carlsberg Meridian Telescope (CMT) has a 2k by 
2k CCD camera with a Sloan r filter operating in a 
drift scan mode. The CMT maps the sky from La Palma 
(Spain) covering the declination range —30° to +50° 
with a magnitude limit of rcMC = 17. The Carls- 
berg Meridian Catalogue, Number 14 ( Version 1.0: CMC 
ICopenhagen University Obs. et aJl I2006T ) is an astrometric 
and photometric catalogue of 95.9 million stars covering 
9 < rcMC < 17. We cross-matched the CMC catalogue with 
2MASS using a 2" matching radius. Because the surveys 
used here avoid the Galactic plane, the contamination from 
matching to other stars within 2" will be relatively small 
l|Girven et al.ll201lh . With these combined catalogues, we 
were able to calculate an (rcMC — J) colour as a diagnos- 
tic for spectral type, as well as (J — K), indicative of strong 
companion star contributions in composite systems (see Fig- 
ure[T}. 

Based upon the (rcMC — J) colour of the composite 
models described in Section^ the CMC sample was cut to 
include only stars bluer than a GOV star (5750 K on the 
ICastelli fc Kuruczl 120031 grid), i.e. (rcMC — J) < 0.9. The 
cut includes all possible combinations of subdwarf plus com- 
panion, but removes a significant fraction of contaminants. 
This does not limit our selection of subdwarfs with compan- 
ions as discussed in Section[4] The sample was also limited 
to rcMC < 16.0, primarily to match the magnitude limit of 
2MASS (K s ~ 14.3). This resulted in ~ 1.9 million objects. 

All objects within the (rcMC — J) colour cut were cross- 
matched with the Gala xy Evolution Explor er (GALEX) all- 
sky ultraviolet survey l|Martin et al.l l2005h Data Release 6. 
This provides magnitudes in two bandpasses, mpuv and 
niNuv, centered around 1500 and 2300A, respectively. 
The matching was performed using the predefined cross - 
matching tables in GALEX CasJobs (jBudavari et al.ll2009l ) 
searching for all sources within 2". The resulting catalogue 
of neighbours contains approximately 560, 000 matched ob- 
jects and hereafter will be referred to as the "C2M" sam- 
ple. The mean of any multiple GALEX observations was 
taken where available an d both bands w e re co rrected for 
non-linearity according to lMorrissev et all (|2007T l. 

Finally, the objects from the match between CMC, 
2MASS and GALEX were further cross-matched with the 
Sloan Digital Sky Su rvey (SDSS) Data Release 7 (DR7; 
lAbazaiian et al.l I2009T ). This sample will hereafter be re- 
ferred to as the U C2MS" sample which is smaller and pho- 
tometrically dee per. The SDSS Ca sJobs predefined cross- 
matching tables (|Li fc Thakarlliooj) were utilised. Objects 



were limited to have good quality photometric magnitudes 
(see Table[T}. This resulted in a sample of ~ 105,000 ob- 
jects for which good SDSS u, g, r, i and z magnitudes were 
available along with GALEX, CMC and 2MASS photome- 
try. For ~ 1.5 percent of objects within this sample, SDSS 
optical spectra are available. 

3.2 Sample II: GALEX, SDSS and UKIDSS 

The GALEX, CMC and 2MASS cross-matched sample dis- 
cussed above benefits from covering a large area (limited by 
the GALEX footprint), but is relatively shallow with a lim- 
iting magnitude of r = 16.0 and Ks — 14.3. This restricts 
our ability to construct volume-limited samples. 

A second, complimentary sample was selected from 
GALEX, SDSS and UKIDSS. One of the five UKIDSS sub- 
surveys, the Large Area Survey (LAS), aims to be the in- 
frared counterpart to the SDSS. UKIDSS LAS will even- 
tually provide imaging over 4028 deg 2 in four broad band 
colours, Y, J, H, and K, with limiting (Vega) magnitudes 
of 20.2, 19.6, 18.8 and 18.2, respectively. This adds a sig- 
nificant increase in depth over 2MASS. Here, w e made use 
of UKIDSS Data Release 9 (see lDve et al]|2006| ). which cov- 
ers ~ 60 per cent of the total, planned, LAS area. SDSS and 
UKIDSS were cross-matched to find the closest match within 
2" using the UKIDSS-SDSS pre-match tables. This sample 
was then matched to GALEX within 2", using the CasJobs 
neighbours search, returning approximately 120, 000 ob- 
jects. Again, multiple GALEX neighbours were combined 
into a sin gle measurement and fl uxes were corrected for non- 
linearity IjMorrissev et al.l 120071 ). This sample will hereafter 
be referred to as the "SU" sample. It is limited in area by 
the current size of UKIDSS, but extends several magnitudes 
deeper than 2MASS in K. Because the UKIDSS LAS area 
is entirely encompassed by the SDSS footprint, we can make 
use of the higher precision, deeper SDSS photometry, rather 
than CMC. The number of objects at each stage of the anal- 
ysis is given in Tabled 



4 SELECTING ULTRAVIOLET EXCESS 
OBJECTS 

4.1 Colour-colour diagrams 

Figure[T] shows colour-colour diagrams for the objects with 
detections in GALEX, CMC and 2MASS. We compare 
(hifuv - rcinc) vs (rcMC - and (m N uv - rcMc) vs 
(rcMC — K a ), where the (J — K s ) colour of each object is 
colour-encoded in the plot. For single stars, the (J — K„) 
range corresponds to spectral types 05 to K0. The colour 
indices are tailored to highlight in colour-colour space the 
position of composite blue plus red objects. The (rcMC—K s ) 
colour of an object is a relatively good indication of stel- 
lar spectral type and (mFuv — rciuc) will indicate objects 
with an excess in the ultraviolet in contrast to single main- 
sequence stars. The truncation at (rcMC — K s ) ~ 1.5 is 
caused by our imposed cut of (rcMC — J) < 0.9. 

In Figure[2] the same sources are plotted but now encod- 
ing the density of sources on a grey scale to better represent 
relative numbers. The main sequence is found along the bot- 
tom edge of the main group of objects in the (mpuv — rciic) 



© 0000 RAS, MNRAS 000, 000-000 



J.Girven et al. 




Figure 1. Colour-colour diagrams of (mpuv — »"cmc) vs (''CMC — K s ) (left; ~ 15,000 objects) and (niNUV — fCHc) vs 
(rcMC — K a ) (right; ~ 105,000 objects) for objects with matches in GALEX, CMC and 2MASS catalogues (the C2M sam- 
ple) and satisfying the magnitude cuts in Tablefl] The majority of objects are main-sequence stars and therefore do not have 
mpuv magnitudes . This causes the fac tor of 10 difference in numbers of objects in the two diagrams. E(B — V), taken from the 
GALEX catalogue Schle gel et al. I Jl998ft map value, is also limited to < 0.15 and for clarity, mnjuv *S 16.5. The plotted colour for 
each point corresponds to the (J — K s ) value, where colours above 0.5 or below -0.3 are shown a s 0.5 and -0.3, respecti vely. For single 
stars, this range corresponds to spectral types 05 to K0. The sdB/sdO star candidates from iKilkennv et al.l l ll988h with matches 
in the C2M catalogue are shown as black circles (totalling 84 objects). A reddening vector corresponding to E(B — V) = 0.15 is 
shown as a black arrow centered on (0, 5) in both diagrams. The dashed black lines show the colour-colour selections as discussed 
in Sectionl4~2land TablcQ] 




Figure 2. 2D density plots of the (mpuv — ^CMc) vs ('"CMC — Kg) and (mmjv — f*CMc) vs (''CMC — Ks) colour-colour diagrams 
in FigureTTI however, hintjv ' s no longer limited to 16.5. The grey scale is shown on the respective colour bars, labelled by the 
number of objects (N) on the left and N/1000 on the right. The reddening vector again corresponds to E(B — V) = 0.15 and 
is shown as a black arrow centered on (0,5). The dashed black lines show the colour-colour selections as discussed in TablefT] 
Additionally, subdwarf (15,000 ^ T e g ^ 40, 000 K) and main-sequence star model colours (Section[2} are plotted as blue and red 
open circles, respectively. The yellow lines show a single subdwarf (15, 000 ^ T e ff 40, 000 K in steps of 5, 000 K) paired with the 
sequence of main— sequence stars described in Section[2] As a reference, subdwarfs (of the above temperatures) with approximately 
a K0 or FO-type companion are joined with a black line. 
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Table 1. Colour selection for finding subdwarfs with companions, for both the C2M and C2MS (with or without SDSS magnitudes) 
and SU samples. Constraints with a "Sample" flag were only applied to that sample, whereas constraints with no flag were applied to 
both samples. "bad_flags" is defined as saturated or bright or edge or nodeblend and "nChild" is the number of children objects detected 
by SDSS. 



Colour 


Constraint 




Sample 


rCMC 




16.0 


C2M 


(rCMC - J) 


< 


0.9 




rpMC Uncertainty 




0.10 




E(B — V) 




0.15 




( niT?T T~\? — TCKKC" I 
\ r U V \^lvl\^/ 




3.8 * (rn-\/m — Ki) — 0.3 




(m FUV — r C MC ) 




-2.7 * (rcMC - K a ) + 4.7 




(mpuv — ''CMC) 




-3.0 




FUV artifact flag 




1 




niFUV Uncertainty 




0.05 


C2M 






0.10 


SU 




AND 






(m NUV - r CM c) 




1-3 * (rcMC - K s ) + 0.54 




(m NUV - rcMc) 




-1.45 * (r C MC - K s ) + 3.3 




(m NU v - fCMc) 




3.5 * (r CMC -K s ) + 0.12 




(m NUV - t"CMc) 




-2.0 




NUV artifact flag 




1 




(rCMC - K s ) 




1.75 


C2M 


niNUV Uncertainty 




0.05 


C2M 


(''SDSS - K) 




1.5 


SU 


m NUV Uncertainty 




0.10 


SU 


SDSS specific: 








flags & bad_flags 









nChild 










Table 2. Summary of numbers at each stage of the processing. The left hand columns shows which surveys were included at that stage 
in the processing. The "In cuts" columns are those objects satisfying the criteria from Figure[T]and Table[T]and the "S" column displays 
the total number of objects in this catagory. 



Sample 


Surveys 






Total 




In cuts 




Name 


Ultraviolet 


Optical 


Infrared 


(approx) 


S 


SDSS Spectra 


SIMBAD 






CMC 


2MASS 


1,900,000 








C2M 


GALEX 


CMC 


2MASS 


560, 000 


449 




58 


C2MS 


GALEX 


SDSS CMC 


2MASS 


105, 000 


93 


25 


21 






SDSS 


UKIDSS 


220, 000 








SU 


GALEX 


SDSS 


UKIDSS 


120, 000 


134 


72 


47 



vs (rcMC — Ks) plane, and more centrally through the main 
group in the (m N uv - rciuc) vs (r C MC - K 3 ) plane (Fig- 
ure[T}. Simulated colours derived from our main-sequence 
star model (Section[2| confirm that this is the expected po- 
sition of the main sequence in our chosen colours. Similarly, 
composite sdB plus companion star models are also shown 
in Figure[5] highlighting the region of colour-colour space 
where we expect to find such systems. 

The large scatter in (m FU v - tcmc) or (m N uv — tcmc) 
for a given (rcivrc — K s ), especially at the red end, can be 
explained due to a few factors. First of all, even though 
we formally require the tcmc uncertainty to be less than 
0.1, there appears to be additional systematic scatter in the 
rcMC magnitudes. Investigating the vertically extended re- 
gions in our colour-colour diagrams (Figure[T]) when using 
the much more reliable rsDSS instead of rcnic, we find that 
the spread is significantly reduced. However, the larger sky 



coverage of the CMC is far more important for our study 
especially as the subdwarf plus companion systems fall in a 
relatively clean part of the diagram. Another reason for the 
observed spread is the fact that the GALEX magnitudes 
have been shown to suffer from non- linearities for bright 
stars, amongst o ther problems (e.g. iMorrissev et al.l 120071 ; 
I Wade et alj|2009h . Although w e corrected for non-line arity 
using the method described in IMorrissev et ail (120071 ). the 
equations are empirical and there may be a significant scat- 
ter in individual measurements. 

In addition, despite limiting E(B — V) ^ 0.15, much 
of the spread around the main sequence can be accounted 
for by considering the effects of interstellar reddening. The 
E(B — V) magnitude for each object is taken from the 
GALEX catalogue, w hich is itself calculate d from the Galac- 
tic reddening maps of lSchlegel et al.l (|l998l ). The interstellar 
reddening is illustrated by the reddening vectors in Figure[T] 
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T hese are calculated by foldin g the mean extinction curve 
of iFitzpatrick fc Massal (|2007f l through the relevant filter 
transmission curves. In the (mpuv — rcinc) vs (rcivic — K s ) 
plane, reddening of blue objects moves them above the main 
sequence in (hifuv — rcMc), a region populated by a num- 
ber of objects. However, reddening in the (mNuv — rcMc) 
vs (rcMC — K s ) plane approximately moves objects along 
the main sequence. The components of the reddening vec- 
tors are approximately the same in both (mFuv-fcinc) and 
(mNuv — rc Mc) because of the 2200A b ump in the redden- 
ing function |Papoular fc Papoular|[2009l ) coincides with the 
central wavelength of niNuv- However, the intrinsic and sig- 
nificant variations in the reddening law along different lines 
of sight affect the ultraviolet magni t udes m ore so than the 
optical values. IFitzpatrick fc Massal 1)20071 ) show that even 
when considering the standard stars that are used to cal- 
culate the adopted reddening function, a significant spread 
around the mean extinction curve is observed. This leads 
to large departures from the mean law, affecting the ul- 
traviolet region in particular. These variations in the ex- 
tinction curve, along with the variation of the true redden- 
ing to the subdwarf compared with that calculated in the 
ISchlegel et all <|l99ct ) maps, are thus likely responsible for 
the stellar sources populating a vertically extended region 
in the (hinuv — rcMc) vs (rcMC — K s ) plane. In any case, 
the outliers form only a small fraction of the total source 
population and the reddening vector does not move main- 
sequence stars into the colour selections we discuss below. 

4.2 Isolating subdwarfs in binaries 

In order to classify our sources and check for known ob- 
jects within our sample, we resolved all positions using 
simbaeQ, and also consulted any available SDSS opti- 
cal spectra. In the upper-left corner of the (hifuv — rcnic) 
vs (rcMC — K s ) colour-colour diagram, one would expect 
to find white dwarfs and single-star subdwarfs, which is 
corroborated by classifications in the SIMBAD database. 
Unfortunately, none of our sources with colours consistent 
with single subdwarfs have SDSS spectra that could conclu- 
sively confirm their classification (due to them saturating in 
SDSS). The objects towards the right of the diagram, with 
(rcMC — K) ~ 2.0, prove to be galaxies. These are removed 
by use of the point s ource flag in SDSS. 

iKilkennv et al.l (1 19881 ) created a catalogue of subdwarf 
stars and candidates from previous studies, including work 
on the PG survey. This includes subdwarfs both with and 
without companions. We matched this catalogue to the C2M 
catalogue, resulting in 1704 objects. The subset for which ap- 
propriate quality limits are satisfied are plotted in Figure[T] 
(84 sources). We see that this sample splits into two dis- 
tinct groups. A significant fraction falls in the region where 
single subdwarfs and white dwarfs are expected to lie. How- 
ever, a good number of these (~ 35 per cent) lie at a much 
redder (rcMC — K s ) colour, where, from the synthetic mag- 
nitudes calculated in Section[2l we expect subdwarfs with 
main-sequence star companions. The objects in this redder 
region (inside the black dashed lines in Figures [T] and [5]) , 
would appear to be main-sequence F or G-type stars from 

1 http: / /simbad. u-strasbg.fr/simbad/ 



their (rcMC — K s ) colour, but have an ultraviolet excess in 
(m FU v - rciic) and/or (m N uv - rc Mc) colour. Th i s con - 
firms that a significant fraction of the IKilkennv et al] (|l988T l 
sample show photometric evidence for being composite, but 
also that we have detected a large number of new sources 
within that same region of colour space. 

For the new C2M objects in this region, where SDSS 
spectra are available, they can be seen to be mostly subd- 
warfs along with one white dwarf and two cataclysmic vari- 
able stars (CV: see Table[3}. SIMBAD, however, only re- 
turns four known subdwarfs in this region of colour-colour 
space. This may be expected as previous work has intention- 
ally focused on single-lined sdB systems that are therefore 
dominated by the subdwarf. The number of objects grouped 
under a few broad classifications are summarised in Table[3] 
Note that close to 90 per cent of the C2M sources (without 
SDSS) within this region are unknown. 

In order to isolate composite subdwarfs while avoid- 
ing obvious contaminants, we devised cuts in colour-colour 
space (Table[T]) guided by our simulated composite sub- 
dwarf colours and the SIMBAD and SDSS spectroscopic 
classifications discussed above. The right hand side of the 
cuts was chosen to avoid contamination from galaxies and 
quasars, and similarly on the lower side the main sequence 
was avoided. At the left hand edge, the cuts were chosen 
to avoid early-type stars and single subdwarfs. We require 
objects to be in both the (mFuv — rciic) vs (rcMC — K s ) 
and (mNuv — rcMc) vs (rcMC — K s ) cuts because objects 
residing in just an individual box are likely to arise from 
spurious GALEX fluxes. Contamination of this region due 
to interstellar reddening is small because very few objects 
will be moved from the main sequence, along the redden- 
ing vector, into the box, as shown in Figure[T] Similarly, the 
scatter from a poor rcMC magnitude does not lead to a large 
contamination, because the subdwarfs with companions re- 
gion is sufficiently far from the main sequence. 

We repeated a similar selection using the SU sam- 
ple, again using (m FU v - t S dss) vs (r S Dss - K) and 
(mNuv — tsdss) vs (rsDSS — K) colour-colour diagrams (not 
shown). All magnitudes were limited to have uncertain- 
ties less than 0.1 mag and E(B — V) ^ 0.15. An increase 
in the number of quasars was seen, which encroached on 
the cuts used for (niNuv — rcinc) vs (rcMC — K s ). The up- 
per limit on (rsDSS — K) was therefore reduced, as shown 
in Table[T] however the contamination was not completely 
removed. The cuts on (mFuv — ?"sdss) vs (rsDSS — K) re- 
mained unchanged, where we ignore the small differences 
between UKIDSS K magnitude versus 2MASS K s magni- 
tudefl After these adjustments, 134 objects reside within 
the cuts, 72 of which have SDSS spectra. This is signifi- 
cantly more than the C2MS sample because many of the 
C2MS objects are saturated in SDSS. As for the 2MASS 
sample, we provide broad classifications for the SU sample 
in Table[3] 

With our selection cuts in place, we can use the tracks 
of our synthetic subdwarf-companion pairs to consider the 
completeness of our composite subdwarf sample. We find 



2 Assu ming a J — K 3 co lour of ~ 0.3 and using the transforma- 
tions of [Carpenter (200l|), the difference between the K 3 and K 
magnitude is ~ 0.003 and therefore negligible. 
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Table 3. Table of classifications for the 449, 93 and 134 objects in the C2MS, C2MS and SU samples, respectively, and inside the 
colour-colour selection boxes from Figuro[T]and Table[T] The SDSS spectra column is from visual inspection of the optical spectra. The 
number of galaxies seen in the SDSS spectra is virtually zero because the flags used to select the SDSS objects remove any extended 
objects. 





C2M 




C2MS 




SU 


Classification 


SIMBAD 


SIMBAD 


SDSS spectra 


SIMBAD 


SDSS spectra 


SD 


7 


4 


22 


7 


62 


Composite 


9 


1 





4 





CV/Nova 


21 


8 


2 


10 


4 


Galaxy 


2 











2 


Quasar 

















WD 


19 


11 


1 


26 


4 


Total with classification 


58 


24 


25 


47 


72 


Total without classification 


391 


69 


68 


87 


62 



that our region covers only a limited range in companion 
type for a given subdwarf temperature, as systems that are 
either dominated by the companion or the subdwarf fall out- 
side our region. This choice is required to reduce contami- 
nation from single stars. Based on our simulated colours, 
we find that subdwarfs with temperatures up to 30, 000 K 
would fall in the (hifuv — rcMc) vs (rcivic — K s ) colour cut 
for even the coolest main-sequence companion in our grid 
(3,000K: ~M5). 35, 000 K and 40, 000 K subdwarfs, how- 
ever, would require > 3, 750 K (<M0) and > 5, 000 K (< K0) 
companions, respectively, to make them stand out from the 
main sequence populations. In the case of early-type com- 
panions, subdwarfs plus O-type and B-type stars are also 
lost as they merge back into the blue end of the main se- 
quence. A 15, 000 K, 20, 000 K, 30, 000 K and 40, 000 K sub- 
dwarf would be identified if it had an < 7, 500 K (>F0), 
< 8,250K (>A5), < 8,000K (>A5) or < 8,500K (>A5) 
companion, respectively. 

For the colour-colour tracks, the companions are re- 
stricted to be main-sequence stars. However, we may also 
expect to find a population of subdwarfs with sub-giant 
or gi ant companio ns similar to HP 185510 (|Fekel fc Simon! 
Il985l). HP 128220 jHowarth fc Heberlll990r i and BP- 7° 5977 
l|Viton et ail Il99ll ; iHeber et all 12002V In fact, the binary 
population synthesis of Han et all (|2003T ) predicted that 
the majority of K-type companions to subdwarfs should be 
evolved companions. We calculated the (mpuv — ?"cmc) vs 
(rcMC — K a ) location of G7 to K3-type giant stars (Figure(3] 
upper-right panel) by taking the solar metalicity, zero age 
horizontal branch stars from the ICastelli fc Kuruca l|2003h 
model atmosphere library and again rescaling the fluxes to 
a corresponding z ero age horizontal b ranch luminosity from 
the isochrones of iGirardi et al.l (|2000l ). All combinations of 
subdwarf plus giant star systems fall outside of the colour 
cuts described in Table(3] Systems with either overluminous 
subdwarfs, or companions in an intermediate state between 
the main-sequence and the horizontal branch may, however, 
fall within the colour cuts. We do not expect these to be a 
significant population in our sample. Since we do not expect 
specific formation mechanisms to become more or less preva- 
lent as a function of distance, we can still use our sample to 
study the spatial distribution of subdwarfs even if the sub- 
sample of subdwarfs with evolved companions is selected 
against. 



We may also expect that some detached white dwarf 
plus main-sequence type companion systems are found to be 
contaminants of the sample, since these are composite sys- 
tems with a hot component and a cooler companion. How- 
ever, we simulated the colours of such systems and, with the 
exception of very low gravity white dwarfs, they do not fall 
in the colour-colour region selected in Table[3] (see Figure(3] 
lower two panels). In this colour space, the small radius of 
the white dwarf means that the flux is dominated by all 
but the latest of main sequence companions and so they lie 
closer to the main sequence in both diagrams. They are thus 
unlikely to constitute a significant contaminant. 

Looking out of the Galactic plane to distances of over 
1 kpc, we may expect to see a sizeable fraction of thick 
disk and halo stars. Therefore, the companions to the sub- 
dwarfs in our samples may be metal-poor. In Figure[3] 
(upper-left panel), we show that subdwa rfs with metal-poor 
(log([M/H]/[M/H] so i ar ) = 1.5: ATLAS9: ICastelli fc Kuruczl 
12003! ) companions indeed still fall in our colour selection. We 
discuss the associated possible biases on our fitting technique 
in Section[6] 

A summary of our sample sizes at various stages of the 
analysis can be found in Tabled The full list of 449 objects 
inside our C2M sample can be found in Tabled 



5 SPECTROSCOPIC OBSERVATIONS 

We discuss here some spectroscopic follow-up obtained to 
verify that the C2M sample objects likely contain a sub- 
dwarf component before turning to the modelling of their 
spectral energy distributions (SEP) in Section [S] 



5.1 WHT 

Nine objects falling within the colour-colour cuts described 
in Table[T]were observed in July and Pecember 2010, using 
the 4.2m William Herschel Telescope (WHT) at the Roque 
de los Muchachos Observatory, La Palma, Spain. We used 
the ISIS dual-beam spectrograph mounted at the Cassegrain 
focus of the telescope, with a R600 grating on both the blue 
and the red arms, and a 1" slit. The blue arm of the spec- 
trograph is equipped with a 2048 x 4096 pixel EEV12 CCP, 
which we binned by factors of 3 (spatial direction) and 2 
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Figure 3. Potential contaminants of the subdwarf plus main— sequence star sample, following the same format as Figure[2] See Sec- 
tion ll.2l for a discussion. Top-left: Subdwarfs (blue open circles: described in Section[2} with metal-poor (log([M/H]/[M/H] so i ar ) = 
1.5: red open circles) main-sequence star companions. Composite objects are shown in yellow. Top-right: Subdwarfs with giant 
star companions (ranging in spectral type from approximately G7 to K3). Single G and K-type stars do not fall in the range of the 
Figure, and therefore we mark their (rcMC — Ks) position by downward pointing red triangles. Bottom-left: log(g) = 8 DA white 
dwarfs (green open circles) with main-sequence companion stars (described in ScctionO- Bottom-right: log(g) = 7, and therefore 
larger radii, DA white dwarfs with main-seq uence companion stars. The DA white dwarf model grid was kindly provided by D. 
Koester (for a description, see|Koestcr 201(3), and ranges from 15,000 ^ T c ff ^ 40,000K in steps of 5,000K. 



(spectral direction). The 2048 x 4096 pixel REDPLUS CCD 
on the red arm was binned similarly. This setup delivers 
a wavelength coverage of 3772 - 5136A on the blue arm, 
with an average dispersion of 0.88 A per binned pixel, and 
5983 — 7417A on the red arm, with an average dispersion of 
0.98A per binned pixel. We determined the resolution to be 
1.2A, from measurements of the full width at half maximum 
of night-sky lines. The setup during the December observa- 
tions was identical, except that the CCDs were binned 2x2. 

The spectra were debiased and flatfielded using the 
STARLlNp|3 packages KAPPA and FI GARO and the n optimally 
extracted using the PAMELA code ()Marsh|[l989l ). We derive 
the wavelength calibration from Copper-Neon and Copper- 
Argon arc lamp exposures taken during the night, selecting 

3 Maintained and developed by the Joint Astronomy Centre and 
available from http://starlink.jach.hawaii.edu/starlink 



the arc lamp exposure nearest in time to each science spec- 
trum. 

Finally, the raw spectra were converted to flux units and 
the telluric absorption lines removed. For the July run, the 
flux calibration was done using a model spectrum of a "flux 
standard" DA white dwarf, observed on the same night. The 
December run suffered from poor weather and no flux stan- 
dard was observed. We calibrated these two spectra using 
an earlier observation of SP1446+259, taken with the same 
instrumental setup. The shape of the spectrum is therefore 
reliable, but the absolute flux level is not. Our analysis does 
not depend on the absolute flux of the targets, so our con- 
clusions are unaffected. 

We plot the resultant spectra in Figure!?] and find that 
all but one of the nine objects chosen from the colour-colour 
selection are sdB stars with companions (Tabled . 
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Table 4. Follow-up spectroscopic observations and classifications. 
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Figure 4. WHT optical spectra of nine candidate subdwarf plus companion stars chosen from the colour-colour selection seen in 
Figurc[T] Eight of the nine targets are subdwarfs with hints of absorption lines from the companion star (see also Section|4]l. Spectra 
are ordered approximately by effective temperature of the subdwarf, and shifted in flux by appropriate amounts. 0051—0955 is 
probably an A-type star rather than a subdwarf. 



5.2 MagE 

In addition to the WHT spectra, two candidates were ob- 
served on 7-8 June 2010, using the MagE (Magellan Echel- 
lette) spectrograph mounted on the Magellan-Clay telescope 
at Las Campanas Observatory, Chile. We used the 1" slit 
with the 175 lines/mm grating to cover ~ 3100 — 11200 A at 



a resolution of R = 4100. The data were unbinned and we 
used the slow readout mode. 

The spectra were reduced with the Carnegie pipeline 
written by D. Kelson. This Python-driven pipeline performs 
typical calibrations: flat-fielding, sky background subtrac- 
tion followed by optimal extraction and wavelength cal- 
ibration. The wavelength calibrations were derived from 
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Thorium- Argon lamp exposures taken during the night, 
which provided ample suitable lines over the entire wave- 
length range. The pipeline selects the closest lamp exposures 
in time to each science spectrum. Raw spectra were then flux 
calibrated using a spectrum of the flux standard Feige 110, 
observed at the end of each night. We find that both objects 
observed have spectra consistent with being sdB stars with 
some evidence for a companion. 

This initial exploration of eleven of our candidates thus 
offers strong evidence that we are primarily selecting com- 
posite subdwarf systems with our colour cuts, with a low 
contamination rate. We discuss contamination of our sam- 
ples further in Sections l7.8l and 17.91 



dwarf effective temperatures by a few thousand Kelvin and 
shifting the distribution of companion types later by a few 
hundred Kelvin. If anything, this accentuates the conclu- 
sions we draw in Section[8] 

A final potential bias to our fitting method is that ap- 
proximately 10 per cent of subdwarfs are evolved and there- 
fore will have lower surface gravities a nd bloated ra dii com- 
pared with their unevolved equivalent (|Heberll2009T y Fitting 
a system with an evolved subdwarf using our subdwarf plus 
main-sequence star model grid (described in Section[2]|, we 
would find that the companion star is cooler and the sub- 
dwarf is hotter than the true temperature. However, this 
situation will most likely result in a high minimum x 2 and 
therefore be flagged as a bad fit. 



6 FITTING COMPOSITE SYSTEMS 

To quantify the likely composition of our subdwarf candi- 
dates, we pursued SED fitting exploiting the broad wave- 
length range of the photometric data that is available. The 
subdwarf star dominates the ultraviolet flux while the main- 
sequence companion clearly dominates in the infrared. This 
permits the decomposition of the SED into two components 
at a common distance. In this section we demonstrate that 
good constraints on both the subdwarf and companion star 
effective temperature can be derived from such fits. The ob- 
served magnitudes were fitted with the grid of subdwarf plus 
main-sequence star magnitudes discussed in Section[2] with 
the additional option of having a subdwarf with no compan- 
ion (shown as MS T e g = K in Table[S] onwards) . This was 
performed by minimising a weighted x 2 whilst varying the 
distance, subdwarf and companion effective temperatures. 
Uncertainties were taken from the one sigma contours in 
the \ 2 surface. This fitting was restricted to the sub-samples 
where SDSS photometry is available, since we require multi- 
band optical photometry in order to decompose the SED. 

Reddening from interstellar dust can potentially have a 
significant effect on the shape of the subdwarf SED, espe- 
cially at short wavelengths. It would therefore primarily af- 
fect the inferred subdwarf effective temperature. The slope 
will be flattened and thus a systematically lower effective 
temperature would be found. Without prior knowledge of 
the reddening to the system, this is not easily corrected 
for. To estimate an upper limit for this effect, we calcu- 
late ^the_j^dd^ning atthe position of the subdwarfs from 
the lSchlegel et al.l l| 19981 ) maps and use these values to first 
deredden the magnitudes. Refitting these values gives a sec- 
ond set of system parameters that will, in general, be over- 
corrected for reddening in comparison to the fits without 
any reddening. The true parameters will lie somewhere in 
between these two limits. 

As shown in Figure[3] subdwarfs with metal-poor com- 
panions fall in the colour cuts defined in Table[3] They are 
not a contaminant, but fitting the metal-poor systems with 
solar metalicity models will lead to biased system param- 
eters. Less absorption in the ultraviolet from metal lines 
means the companions will contribute a fairly significant 
amount of flux at short wavelengths. To test the effect of 
this, we fitted the C2MS and SU samples with a grid of sub- 
dwarfs plus meta l-poor (log ( [M /H] / [ M/H] a o i ar ) = 1.5) com- 
panions from the lCastelli fc Kurucj (120031 ) ATLAS9 model 
atmosphere library. This has the effect of reducing all sub- 



7 FIT RESULTS AND INDIVIDUAL OBJECTS 

All the fit parameters for the 93 objects from the C2MS 
sample (Tabled are given in Table flOl Similarly, the 134 
SU objects are shown in Table lTTI We adopt a somewhat 
unusual notation for the upper and lower uncertainties, de- 
noted by the "{" symbol, because the subdwarf and com- 
panion effective temperature uncertainties are strongly cor- 
related. "{" indicates the upper and lower la uncertain- 
ties added to the best fit value. The upper values all cor- 
respond to the same fit solution and similarly for the lower 
values. As an example, consider a hypothetical system where 
SDT cfi = 15,000{™K, and MS T cff = 2, 000{J;°°° K. 
This corresponds to three solutions: the best fit (a 15, 000 K 
subdwarf with a 2,000K companion), a la uncertainty in 
the direction of increased subdwarf temperature (a 25, 000 K 
subdwarf with a 1,000K companion), and a la uncer- 
tainty in the direction of decreased subdwarf temperature (a 
10, 000 K subdwarf with a 3, 000 K companion). One cannot 
mix and match these combinations. For example, a 10, 000 K 
subdwarf with a 1, 000 K companion, or a 25, 000 K subdwarf 
with a 3, 000 K companion, are not valid solutions. A mini- 
mum uncertainty is set at one grid point and therefore is also 
limited by the extent of the grid: a minimum and maximum 
subdwarf temperature of 11, 000 and 40, 000 K, respectively. 
We examine systematic uncertainties in Section ing leading 
to estimates of a few thousand Kelvin for a more realistic 
error. We show example SEDs and fits to a few objects in 
Figures[S]and|Sl Objects in Figures[S]and|§]are found to have 
approximately GO and A7-type companions, respectively. 

We compared our results to published effective tem- 
peratures and/or known companions for the C2MS and 
SU samples, shown in Table[5] and [6] respectively. The best 
fit is not always satisfactory, indicated by a high % 2 - We 
include the "Q" (Quality) column to show where this is 
the case. "Q" values correspond to; l:Good fit, 2:Average 
fit, 3:Poor fit, 4: WD/ WD+MS/CV and 5:Quasar/Galaxy. 
Values of three and above are excluded from the histograms 
shown in Figure[8] [9] and [10] The classifications in this 
catagory between values of 1, 2 and 3 are purely qual- 
itative. SIMBAD has an entry for many more objects, 
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2001 
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|2002|, Edelmann et al.1 


2003. 


Morales-Rueda et al.l 


20031. Stark & Wade 


2003. 


Naoiwotzki et alj|2004 Reed & Stienind 12004, Lisker et al. 


20051, 0stenserJl2OO6l,|Wade et al.|l2006l.lStark & Wade 


2006. 


Stroeer et all 120071, Iwade et all 120091, iGeier et all 1 


2011a 


and Vennes et al. 12011 


to be composite subdwarf plus 



7.3 1300+0057 and 1538+0934 

The published effective te mperatures f or 1300+0057 
(39359 K: HE 1258+0113: IStroeer et al 



2007) 



companion systems are highlighted in TablesfTUl and [Til 



7.1 Potential systematic temperature differences 

When comparing the system parameters calculated herein 
and those from the literature, there are a number of possi- 
ble causes for discrepancies: Firstly, one must consider the 
fact that often in the literature fitting is performed on the 
absorption line profiles of th e subdwarf with a single star 
model (e.g. ISaffer et al] [l994), whereas our study suggests 
that these systems all have a significant contribution from 
the companion. The single subdwarf fit would then result in 
biased system parameters. 

Secondly, if the subdwarf 's companion is a sub-giant or 
giant type star, our method would underestimate the sub- 
dwarf's effective temperature because we only use main- 
sequence star models for the companion. While this may 
affect isolated cases, we do not expect a significant popula- 
tion of sub-giant and giant companion stars to be present in 
our sample given the colour selection cuts we employed (see 
SectionEOi 

Finally, the suppression of the subdwarf's ultraviolet 
flux due to line blanketing could cause a biased effective 
temperature. Subluminous B stars show peculiar abundance 
patterns. Some metals (mostly the lighter ones) are found to 
be strongl y depleted, while heavier elements can be strongl y 
enriched jo'Toole fc Hebeij 120061 : iBlanchette et al.1 l2008h . 
The abundance patterns are caused by at omic diffusion, 
which depends on various parameters (see iMichaud et al.l 
l201ll . for the state-of-the-art of modelling), however, metal- 
icity may not be an important one. Because the abundance 
pattern differs from star to star, the ultraviolet line blocking 
for any individual subdwarf will deviate from that predicted 
from the solar metalicity models adopted here. Therefore, we 
cannot quantify the systematic uncertaint y in the tempera- 
ture d etermination of the subdwarf stars. lO'Toole fc Heber l 



(2006) regard solar metalicity models as appropriate for sdB 
stars cooler than about 30000 K, but prefer models of scaled 
supersolar abundances for hotter stars as a proxy for en- 
hanced ultraviolet line blocking. Because the effective tem- 
peratures of our program stars are mostly below 30000 K, 
we stay with solar metalicity model spectra. 



1538+0934 (35114 K: HS 1536+0944: iLisker et al 



an d 
1H), 

both in the SU sam ple, are only uppe r limits on the 
effective temperatures. ILisker et alj ((2005) note the pres- 
ence of a cool (~K0-type) companion in the spectrum 
of 1538+0934 and therefore specifically st ate that the 
estima ted temperature is an upper limit. IStroeer et al.1 
l|2007l) also note the presence of a cool companion based 
on the B — J colour for 1300+0057 and therefore one 
may assume the temperature is also an overestimate. In 
both cases, the best fit model (30,000 and 23, 000 K for 
1300+0057 and 1538+0934, respectively) corresponds to 
a bluer (mpuv — hinuv) colour than the GALEX fluxes. 
Therefore using the higher published effective temperature 
model would not agree with the data. 



7.4 1517+0310 and 1518+0410 

In the case of 1517+0310 and 1518+0410 (PG 1514+034 
and PG 1515+044, respectively: SU sample), the compan- 
ion effective temperatures measured (6, 000{g'ygQ K and 
5, 500{g'25o K, respectivel y) are sign i ficant ly different from 
that in the catalogue of l0stensenl (|2006h (K2 and K4.5; 
corresponding to effective temperatures of ~ 4, 800 and 
4,300K, respectively). The whole SED of 1517+0310 is not 
particularly well fit by the calculated best model. The sys- 
tem has a very blue (hifuv — hinuv) colour and therefore 
the best fit model is forced to be a hot subdwarf, which leads 
to a correspondingly increased companion effective temper- 
ature. 



7.5 1709+4054 

1709+4054 ( PG 17 08+409: C2MS sample), was classified by 
ISaffer et al.1 (|l994l ) to be a subdwarf with an effective tem- 
perature of 28, 500 K. We determined 25, 000 - 28, 000 K if 
we apply no r eddening and 27, 000 — 29, 000 K when ap- 
plying the f ull |Schlegel et al.l (1 19981 ) reddening. However, 
Saff er et al.1 (|l994r 7 fit the line profiles of this composite sys- 
tem with a single star subdwarf model, and therefore com- 
paring the two sets of temperatures is not comparing like 
for like. 



7.2 0018+0101 

ILisker et all {2005) calculated an effective temperature for 
0018+0101 (HE 0016+0044) of 28, 264 ± 800 K. This com- 
pares relatively well with our SU sample estimate of 
25, 000 — 23, 000 K, however, a significantly higher temper- 
ature is measured when using the C2MS sample (39, 000 — 
40, 000 K). Either a 23, 000 or a 40, 000 K subdwarf provide 
an adequate fit to the SED, and small changes in the \ 2 sur- 
face lead to the alternate solution. The flat \ 2 surface comes 
about from a very blue (mFuv — mNuv) colour (—0.64) that 
is difficult to reconcile with the rest of the SED. 



7.6 2138+0442 

For the case of 2138+0442 (PG 2135+045; C2MS 
sample), we find a slightly lower effective temper- 
ature (24,000 - 26, 000 K) compared w ith the pub- 
fished value of lAznar Cuadrado fc Jefferv (~ 28, 000 K: 
I2002D . Including the full ISchlegel et all i| 19981 ) reddening 
(26,000 — 28, 000 K), how e ver, the temperatures agree. 
lAznar Cuadrado fc Jeffervl (|2002l ) treat 2138+0442 as a 
composite system fitting both objects in the blue region of 
the spectrum, thus the above mentioned problem of fitting 
a single star model (Section l7.1|l does not apply. 
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Figure 5. The SEDs of, and fits to, 0316+0042 (PG0313+005), 0814+2019 and 1212+4240 (PG 1210+429). The optical SDSS 
spectra are shown in grey. The GALEX, SDSS, CMC and 2MASS magnitudes are plotted in blue, green, yellow and red, respectively 
with corresponding errorbars. The fit to 0316+0042 comprises a T e fj- = 28, 000 K sdB model and a 6, 250 K star (black dashed 
lines). Similarly, a T cff = 21, 000 K sdB model and a 5, 500 K star is used for 0814+2019 and a T e g = 23, 000 K sdB model and a 
5,750K star for 1212+4240. The composite spectra and magnitudes are the solid black line and open black circles, respectively. 
The absolute flux level of the SDSS spectrum does not match the model well in all cases. This is most likely related to calibration 
issues of the SDSS spectroscopy, as it disagrees with the SDSS photometry. For 1212+4240, an archive IUE ultraviolet spectrum 
plotted in grey. 0316+0042 and 1212+4240 are offset in flux by 0.30 and 1.05 units, respectively, for clarity. 



7.7 2244+0106 

2244+0106 (PB 5 146) was foun d to be a post-EHB st ar wit h 
a high velocity in lTillich et all (Hyper-MUCHFUSS; l2CJlll) . 
They estimate a T eff = 33580±680 K, log g = 4.75+0.20 and 
a distance of 18.29 ± 2.45 kpc, compared with our 22, 000 — 
26, 000 K at 6.1 — 7. 9 kpc. However, the companion star is 
not accounted for in iTillich et all (|201ll ) and therefore the 
subdwarfs effective temperature is probably overestimated. 
This is also consistent with the unusually low surface gravity. 



7.8 Overlap 

Where the C2MS and SU samples overlap, a comparison of 
the fits is given in Table[7] and shown in Figure[7| The two 
sets of fits appear consistent within the uncertainties. We 
analysed the distribution of the difference between all the 
C2MS and SU parameters (distance, subdwarf and com- 
panion temperature) and find that the distributions are all 
approximately Gaussian, centered about zero. We do not 
find any evidence to suggest that the two samples effec- 



tive temperatures are systematically offset. The errors on 
the subdwarf effective temperature from the \ 2 fit ma y be 
slightly underestimated, and a more realistic error is a few 
thousand Kelvin. The one difference is that the UKIDSS 
data should better constrain the companion star effective 
temperature due to the greater depth and higher photomet- 
ric accuracy of the near-infrared data. 

Overall, in individual cases, we must bear in mind that 
we may occasionally select the wrong solution (in cases 
where the x 2 surface is relatively flat), nor can we iden- 
tify the exact amount of reddening that should be corrected 
for. However, this study is aimed at providing a statistical 
analysis of the sample rather than correct parameters for all 
individual systems. The errors in the measured parameters 
should be randomly distributed and therefore not effect the 
distributions. It is thus not a significant issue for the analysis 
presented here, but these uncertainties should be considered 
when consulting the fitted parameters of individual systems. 

We saw earlier that the key contaminants in our 
colour box are composite systems containing white dwarfs 
(Table©. Indeed, from the C2MS sample 0018+0101, 
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Wavelength [A] 

Figure 6. The SEDs of, and fits to, 0818-0701, 0825+1202 and 1530+1204, following the same format as Figure^ The fit to 
0818-0701 comprises a T cff = 22, 000 K sdB model and a 7, 750 K star (black dashed lines). Similarly, a T cff = 22, 000 K sdB 
model and a 8, 250 K star is used for 0825+1202 and a T cff = 11, 000 K sdB model and a 8,000K star for 1530+1204. 1530+1204 
and 0825+1202 are offset in flux by 0.40 and 0.85 units, respectively, for clarity. 



0141+0614, 0923+0652, 2117-0015 and 2117-0006 are can- 
didates for being DA white dwarfs with infrar ed excesses 
based on their photometry (|Girven et al.|[201ll ). However, 
such a classification can only be confirmed through follow-up 
spec troscopy. SPSS spec troscopy is available for 2117—0006, 
and iGirven et all |201ll ) classify it as a "Narrow Line Hot 
Star" (NLHS), which they believe to be a group primarily 
made up of subdwarfs. 0018+010 1, discussed in Sectio n lT2l 
is also catalogued as a NLHS bv lGirven et~afl (|201lf ), cor- 
roborating the subdwarf label. 



inants, we know that 12 per cent (3/25) of the objects with 
SDSS spectra are contaminants. Therefore, approximately 
eleven (12 per cent of 93) of the whole C2MS sample will be 
contaminants. The two CVs and one white dwarf with SDSS 
spectra (Tabled can be removed from the histograms. Thus, 
the contamination of the C2MS sample (now with and with- 
out SDSS spectra) used for calculating distributions will be 
9 per cent (8/90). Since any such contaminants will be dis- 
tributed right across our fit parameters, we believe they do 
not distort our statistical analysis to a significant degree. 



7.9 Distributions of fits - C2MS sample 

The distribution of subdwarf and companion effective tem- 
peratures for the C2MS sample is shown in Figure[H]and the 
distribution of distances for the C2MS and SU samples in 
Figure[9] Here we compare the parameters with and without 
reddening corrections. Objects that are known to be contam- 
inants, such as white dwarfs and CVs, have been removed 
from all three (distance, subdwarf and companion temper- 
ature) histograms. For galaxies, these should be flagged by 
SDSS and are therefore removed by the flags in Table[T] 
Using Table[3]to estimate the remaining fraction of contam- 



7. 9. 1 Subdwarf temperature distribution 

Taking the system parameters calculated without correcting 
for reddening, we find that the subdwarf effective tempera- 
tures (Figure[8]) are spread from 20, 000-30, 000 K and peak 
in the 20, 000—24, 000 K range. We do see a pronounced drop 
in numbers belo w 20, 000 K. Reddenin g is not the issue here; 
applying the full Schlcacl et al.l (|l998T l reddening correction 
to the objects before fitting does not lead to a significant 
shift in the distribution, though it is slightly smoothed. 

Based on our theoretical tracks for composite systems, 
we know that we have a reduced completeness below ~ 
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Table 5. Individual objects of interest from the C2MS sample. Only subdwarfs with measured system parameters are displayed, but 
all known CVs and white dwarfs are shown. The "{" notation is described in Section[7]and does not simply represent uncertainties. 
T he fit parameters sho wn are not corrected for interstellar reddening. The comments quoted for the possible white dwarfs matched 
in lGirven et al.l ll201ll) are the classifications according to the SDSS spectra found therein. NLHS corresponds to Narrow Line Hot 
Star (probable subdwarf). 
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sdB T cff 


MS Type Ref / Comments 






(1000K) 


(1000K) 


(kpc) 


(K) 





Subdwarfs 



0018+0101 


HE 0016+0044 


40{| 


5-50{|;J| 


i-5{l:I 


28264 




Lisker et al. (20051 
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Figure 7. A comparison of fits using the C2MS sample versus that using the SU sample, where there is overlap. Objects with a 
"Q"^ 3 in Table[7]are excluded. The dashed line shows a one to one relation between the parameters. 



25, 000 K (see Section l4~2)l . For example, cool subdwarfs of 
~ 15, 000 K with an M-type companion will be missed by 
the colour selection. This could lead to a bias towards hotter 
subdwarfs, which we select over a wider range of companion 
types. In addition, the redder (rcMC~ K s ) colours (Figure[5J) 



means that at the K s magnitude limit (14.3), systems will be 
detected down to fainter rcMC magnitudes. This is however 
offset by the increasing intrinsic brightness of cooler subd- 
warfs (because of decreasing log g and increasing radius). 
We previously discussed a bias towards cooler subdwarfs if 
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Table 6. Individual objects of interest from the SU sample. The "{" notation is d escribed in Section[7] and docs not simply 
represent uncertainties. The comments quoted for the possible white dwarfs matched in lGirven et al.l l l201i] l are the classifications 
according to the SDSS spectra found therein. DA and NLHS correspond to DA white dwarf and Narrow Line Hot Star (probable 
subdwarf), respectively. 
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a Noted presence of a cool companion, therefore temperature is an upper limit, see Section ^. 31 
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Table 7. Comparison of fits using the C2MS sample against that using the SU sample where there is overlap. The "{" notation 
is described in Soctionl7l and does not simply represent uncertainties. The "Q" (Quality) column values correspond to; l:Good fit, 
2:Average fit, 3:Poor fit, 4:WD/WD+MS/CV and 5:Quasar/Galaxy. 
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6.50{|:i 


l-5(?i 


1 


29(i 


6 50{|;I| 
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our assumption of a main-sequence type companion is incor- 
rect (Section l7.1[) . However, we believe this to be a relatively 
small fraction given our sample selection (Section l4.2|) . 

To quantify these possible biases, the limitations on dis- 
tance introduced by various magnitude cuts can be seen in 
Table[5] These are derived by taking the absolute magni- 
tudes of the composite system and calculating the distance 
the object would have to be moved to in order to have an 
apparent magnitude at the relevant limit. The primary ef- 
fects in this case are caused by the saturation limit of SDSS 
(?"sdss > 14.1), corresponding to a minimum distance, and 
the faint K s magnitude limit of 2MASS (K s < 14.3), setting 
a maximum distance. These significantly depend on compan- 
ion spectral type (see below) and, to a lesser extent, on sub- 
dwarf effective temperature. It can be seen that the imposed 
roMC magnitude limit does not have an effect because the 
K„ limit is always more restrictive. In essence, in the C2MS 
sample, the 2MASS depth limits the volume over which we 
are reasonably complete. 

7.9.2 Companion type distribution 

As discussed in Section l4~2l the way in which we select subd- 
warfs with companions introduces a bias in companion type. 
We expect our selection to be complete for subdwarfs with 
20, 000 ^ T e ff ^ 35, 000 K and companions in the range A5 
to M5-type. Similarly, including the more extreme subdwarf 
temperatures (15,000 ^ T e ff ^ 40, 000 K), we are complete 
for F0 to KO-type companions. The companion type range is 
smaller in the latter case because, for example, a 40, 000 K 
subdwarf with a M5-type companion does not fall in our 



colour selection, whereas a 35, 000 K subdwarf with a M5- 
type companion does. 

The distribution of companion effective temperature in 
Figure[8] ramps up from early spectral types towards ~G0, 
as might be expected from the initial mass function (IMF). 
On the other hand, the subsequent turn over and drop to- 
wards mid-K-type may be a product of our selection bi- 
ases. A 15, 000 K subdwarf with a MO-type companion sat- 
urates in SDSS at d ^ 1.7 kpc and is too faint for 2MASS 
at d ^ 1.7 kpc (Table[8]). Therefore we are not sensitive to 
all subdwarfs with MO-type companions. The best way to 
reduce such biases and test our completion is by probing to 
fainter _R" s -band magnitudes. This was the key motivation 
behind our second sample, using SU which extends several 
magnitudes deeper and reaches K ~ 17.8, though at the 
expense of limited sky coverage. 

7.9.3 Distance distribution 

The calculated distance distribution seen in Figure[5] shows 
a rapid increase towards ~ 2kpc, followed by an extended 
tail. As we discussed previously, the limitations on distance 
due to our magnitude cuts and limits are important and are 
a complex function of subdwarf effective temperature and 
companion type (Table|8|. There are no clean regions where 
all temperatures and companion types are sampled evenly 
to give a complete, volume-limited sample. If one assumes 
that all subdwarfs (independent of temperature and com- 
panion type) are drawn from same parent distance distribu- 
tion, and we select each subdwarf-companion system with 
equal probability, the distribution shown in Figure[9] would 
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Table 8. Limitations on the distance of subdwarf plus main-sequence star candidates caused by the relative magnitude cuts. This 
is calculated for 15,000, 20,000, 30,000 and 40, 000 K subdwarfs and companions with effective temperatures of 3,250K (M0), 
5,000K (K0), 7, 250 K (F0). The important limits considered are; the saturation of rgoss a t 14.1 (therefore a minimum distance), 
the cut made on r^MC a t 16.0 (therefore a maximum distance), the X s -band magnitude limit of 2MASS at 14.3 and the i^G-band 
magnitude limit of UKIDSS at 17.8. 
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represent the true distance distribution. Therefore we would 
be relatively confident that the peak appears at 1.5 — 2.0kpc. 
However, Table[8] does show that for some combinations of 
subdwarf and companion temperature we are no longer com- 
plete at this peak distance. Here again the deeper SU sample 
can provide us with a more complete sample. 



7.10 Distribution of fits — SU sample 

The corresponding distributions of the subdwarf and com- 
panion effective temperature for the SU sample are shown in 
Figure[T0]with the distribution of distances in Figure[9] The 
subdwarf effective temperature distribution is broadly con- 
sistent with that of the C2MS sample, with most subdwarf 
temperatures between 20, 000 — 30, 000 K . It is also similar to 
that shown for uncontaminated sdBs by I Green et al.l (2006, 
Figure 1). To establish a volume-limited sample, we again 
refer to Table[8] where we contrast the impact of the 2MASS 
versus UKIDSS 7f-band limits. The distances sampled are 
significantly larger, though as before dependent on subdwarf 
and companion temperature. Overall, the SU sample should 
be less biased against finding lower temperature subdwarfs 
compared to the C2MS sample (see the second example 
given in Section [7.9.2l a 15, 000 K subdwarf with a MO-type 
companion would now be detected to 8.5 kpc). This does 
not appear to have increased the numbers of low tempera- 
ture subdwarfs found and thus it appears that their absence 
is not due to our sample biases, but represents an intrinsic 
deficit of cool subdwarfs within the subdwarf population. 
Accounting for reddening (as seen in the grey histogram) 
does not have a large effect, although it shifts the calcu- 
lated subdwarf effective temperatures systematically higher 
by 1,000 - 2,000K. 

Comparing the distribution of companion effective tem- 
peratures to the C2MS sample, the SU sample has a larger 
number of objects with early M-type companions. Hence, 
the SU sample overcomes the main limitation found within 
the C2MS sample, the shallow i\s-band data. The increased 
depth of UKIDSS allows us to probe significantly more sys- 
tems with M-type companions, however we still see a deficit 



compared with K-type companions and earlier. This also 
seems obvious from the lack of systems populating the subd- 
warf plus M-type companion region of the colour-colour dia- 
gram in Figure[T] Selecting subdwarfs with companions later 
than ~M5-type is still limited by the colour selection method 
as discussed previously (Section l4.2p . Probing deeper in the 
A'-band does not help for companion types later than ~M5. 
Accounting for reddening has a complementary effect to that 
on effective temperature. As the subdwarfs become hotter, 
the required companion also shifts to higher temperatures. 

We searched for a correlation between subdwarf and 
companion effective temperatures, but none was found at 
a level above the parameter uncertainties. Better statistics, 
from larger samples, are needed to investigate the subtleties 
of population. 

Overall, when considering confirmed subdwarf systems, 
we believe that the fitting method is producing tempera- 
tures accurate to within a few thousand Kelvin and com- 
panion temperatures to within several hundred Kelvin (a 
few spectral types). There is some disagreement between in- 
dividual fit results when compared with the literature. How- 
ever, our principal goal is not to achieve superior parame- 
ters for individual systems. Indeed, more data are required 
to accurately establish parameters for individual systems. 
Our method does appear to be efficient in finding composite 
subdwarf binaries, while our SED fitting is accurate enough 
to allow us to consider the broad statistical parameter dis- 
tributions within our samples. There will be some influence 
from contaminants. However, the numbers of contaminants 
are a relatively small fraction (Table[3]) and wherever possi- 
ble they have been removed from the distributions. 



7.11 A volume-limited sample 

The advantage of using the SU sample is that significantly 
larger distances are probed. Referring to Table[8] the sample 
is complete for F0 to MO-type companions over distances of 
2.4 to 3.1kpc. Although we can therefore construct a volume 
limited sample in this region, only 14 objects with good qual- 
ity fits to the SED fall within this region (Figure[l0}. It is 
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Figure 8. Distributions of the subdwarf (left) and companion (right) effective temperatures calculated from the fitting method 
described in Section|6] when applied to the C2MS sample. The grey and black histograms show the system parameters when 
calculated with and without the (maximum) reddening correction, respectively. Adjoining pairs of histogram show the number of 
objects in the same bin. A total of 66 objects are included in the histograms, where 27 objects that are known to be contaminants 
(from their SIMBAD classification or their SDSS spectra), or the subdwarf-companion model provides a bad fit ("Q"SS 3 in 
Tablc fTOt . have been removed. The subdwarf effective temperature histogram is grouped in bins of 2,000K and the companion 
star histogram uses bins of 500 K. 
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Figure 9. Distribution of the distance to the subdwarf-companion star systems as calculated from the fitting method described 
in SectionlBlwhen applied to the C2MS (left) and SU samples (right). The grey and black histograms show the system parameters 
when calculated with and without the reddening correction, respectively. Adjoining pairs of histogram show the number of objects 
in the same bin. Objects that are known to be contaminants have been removed. A total of 66 objects are included in the C2MS 
histogram (left), where 27 objects that are known to be contaminants (from their SIMBAD classification or their SDSS spectra), 
or the subdwarf-companion model provides a bad fit ("Q"^ 3 in Table lTOl or II 11 1 . have been removed. In the SU histogram, 84 
objects are included, where 50 have been removed. The distances are given in kpc and the bin sizes are 0.2 and 0.5 kpc for the 
C2MS and SU samples, respectively. The vertical lines in the right hand plot show the region where the volume limit sample is 
defined (2.4 to 3.1 kpc). 
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Subdwarf T eff [ 1 OOOK] Companion [ 1 OOOK] 

Figure 10. Distributions of the subdwarf (left) and companion (right) effective temperatures calculated from the fitting method 
described in Section|6]when applied to the SU sample. The grey and black histograms show the system parameters when calculated 
with and without the reddening correction, respectively. Adjoining pairs of histogram bars show the number of objects in the same 
bin. A total of 84 objects are included in the histograms, where 50 objects that are known to be contaminants (from their 
SIMBAD classification or their SDSS spectra), or the subdwarf-companion model provides a bad fit ("Q"> 3 in Table HCTt . have 
been removed. Similarly, the green and yellow histograms shows distributions (with and without reddening corrections) when the 
distances to the objects are limited to be between 2.4 and 3.1kpc, such that the histogram is a volume limited sample. 11 and 16 
objects are included in the histograms, respectively. The subdwarf effective temperature histogram is grouped in bins of 2, OOOK 
and the companion star histogram uses bins of 500 K. 



impossible to draw meaningful conclusions about parameter 
distributions for such a small sample. 

Following the assumptions described in Section l7,9.3l 
and thus assuming the distribution seen in Figure flOl is rep- 
resentative of the true distance distribution, the peak at 
2 ± lkpc in our distance distribution may then be associ- 
ated with representing the spatial distribution of the bulk 
of the subdwarf binaries. Then, assuming the distribution 
follows a simple disk population of the form r 2 exp(—r/H), 
where r is the distance from the center of the disk and H is 
the scale height, the turning point in a distance histogram 
should represent 2H. Therefore, the scale height of the sub- 
dwarf population in the SU sample is 1 ± 0.5kpc. 



8 DISCUSSION 

Existing samples of subdwarfs have s hown that a subst antial 
fraction of them reside in binaries. Irian et all ([2003) used 
population synthesis models to calculate that the intrinsic 
binary fraction should be 76 — 89 per cent. Our samples ex- 
plicitly target composite systems and thus should be domi- 
nated by subdwarfs with bound binary companions. iHeberl 
(2009) states that the vast majority of subdwarfs have a tem- 
perature between 20, 000 — 40, 000 K. The temperature dis- 
tribution found here appears approximately consistent with 
this range, however we do find a sub-sample of cooler subd- 
warfs with temperatures below 20, 000 K. In the SU sample, 
where sample biases against cooler subdwarfs are smallest, 
they make up ~ 5 — 10 per cent. This is true whether or not 



we account for the full lSchlegel et al.l (1 19981 ) reddening value, 
and thus cannot be an artifact due to reddening. Fitting of 
the ultraviolet part of the SED is especially important for 
calculating reliable subdwarf effective temperature, because 
this is the region where the subdwarf dominates. 

Utilising the SED from the ultraviolet down to the in- 
frared, we have a large range over which both the subdwarf 
and the companion can dominate a region of the spectrum. 
We show that both samples here are sensitive to compan- 
ions of spectral type A5 to M5 for 20, 000 to 35, OOOK sub- 
dwarf effective temperatures and F0 to KO-type if 15, 000 
and 40, 000 K subdwarfs are included. These ranges can be 
seen visually in Figure[2] Many subdwarfs are found to in- 
deed have companions in this regime. In the C2MS sample 
(FigureJH]), the distribution of companion type is seen to be 
a broad peak from F-type companions to ~ KO-type. A sig- 
nificant turnover is then seen towards late K and M-type 
companions, this can be explained, for the C2MS sample, 
because we are only sensitive to these systems over a very 
small distance range. However, the SU sample extends sev- 
eral magnitudes deeper in the 7f-band and therefore removes 
this bias, but still shows a clear deficit of early M-type com- 
panions. This is contrary to the relative abundance of late 
type companions found in many previous surveys. If the M- 
type companions were only paired with cool subdwarfs, they 
would not have been selected by the colour cuts, but this is 
not consistent with the results of the radial velocity stud- 
ies. It therefore appears that subdwarfs with F, G and K- 
type main-sequence companions are intrinsically much more 
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common than those with lower mass M-type main-sequence 
companions, for a broad range of subdwarf temperatures 
(subject to the colour selections described in Sectionl4.2p . 

The population synthesis models of lHan et alj ( |2003l ) 
predict that a significant fraction of subdwarfs will form 
through a channel involving stable Roche lobe overflow. 
These are expected to be ~ 20, 000 K subdwarfs with ~F0 
or KO-type c ompanions close to the main sequence (see Fig- 
ure 15 & 19 o f lHan et alj|2003h . It is believed that these have 
not b een found previo usly because of the "GK selection ef- 
fect" (|Han et alj|2003h . where subdwarfs with F, G and K- 
type companions were not targeted by the PG survey be- 
cause they would show composite spectra (feat ures such as 
the C a II K line and t he G -band). However, IWade et al.l 
\200d\ and lWade et alj l|200gh find that only - 3 per cent of 
the rejected PG stars show indications of being a subdwarf 
with a companion. The majority are (single) metal-poor F 
stars. Here we are primarily selecting subdwarfs with F to 
K-type companions and therefore we would be sensitive to 
this peak. We do indeed find a significant fraction of subd- 
warfs with effective temperatures around 20, 000 K and some 
objects below 20, 000 K. We do not see t he RLOF s ystem s 
dominate quite as strongly as they do in lHan et al.l (2003). 
However, too many cool subdwarfs are found here to be ap- 
propriate for creation solely through the first common enve- 
lope ejection channel (peaking at ~ 30, 000 K) and thus the 
RLOF channel appears to be a significant contributor. 

Part of the iLisker et al] (120051 ) SPY survey sample 
looked at objects with composite spectra. They do not 
find a clear contribution from cool subdwarfs. The SPY 
survey does, however, suffer from strong pre-selection bi- 
ases. The majorit y of targets were sele cted from the Ham- 
burg/ESO survey IjFriedrich et al.l l2000) and required not to 
show evidence of a companion in the low resolution prism 
spectroscopy. The companion types that we are finding in 
this study also a ppear to broadly match the predictions of 
lHan et all l|2003h . The first stable RLOF channel is very 
efficient at producing F to K-type comp anions. Between F 
and K-type companions. lHan et al.l (|2003h predict ~F0-type 
companions to be the most prevalent (by a factor of ~ 3), 
followed by very few ~G0-type companions, and then two 
smaller peaks of approximat ely equal amplitu de at ~ K0 and 
~M0-type (see Figure 15 of lHan et all 120031 ). Our distribu- 
tion does not show the feature at F0, but we may not be 
sensitive enough to FO-type companions, especially in com- 
posite systems with low temperature subdwarfs. Our colour 
cuts only select 15, 000 K subdwarfs with F0 or later type 
companions and therefore we may not show the main peak 
at ~F0-type (right hand panels of FigureslSl and ITT7|) . if it 
is indeed there. Equall y, most of the K- type companions to 
subdwarfs predicted bv lHan et al] (2003) are evolved and lu- 
minous. Therefore they would not be selected in our colour 
cuts because the luminosity of the companion would domi- 
nate the subdwarf. If any of these objects are selected, they 
will be fitted photometrically as a much hotter companion 
than K-type, therefore enhancing the FO-type peak or broad- 
ening it. Thus we should not detect the peak at K0. 

In a simulation, we took a theoretical sample of subd- 
warfs with comp anions that matched the distributions from 
lHan et all (|2003l ) and used the surface gravities discussed in 
Section[6] We then applied the magnitude and colour cuts 
relevant for the C2MS and SU samples. The objects which 



satisfy these criteria do show a similar distribution in effec- 
tive temperature and companion type to that seen in the real 
samples, again suggesting that our observed sam ples are in 
broad agreement with the model populations of lHan et al.l 
(|2003h . 

More recently, IClausen et alj (2012) present indepen- 
dent population synthesis calculations of subdwarfs. In their 
Figure 13, the distribution of companion effective tempera- 
ture is shown using a variety of input model para meters. Run 
6 is th e most comparable to the distribution from lHan et al.l 
( 2003) in terms of input parameters. In this run, and the ma- 
jority of others. [Clausen et al.l (I2012T ) predict a vast majority 
of M-type or later companions to the subdwarfs. This does 
not agree with our samples, which show a lack of M-type 
companions and a significant proportion of K-types. This 
suggests that observational samples such as those presented 
here have the ability to directly constrain binary population 
synthesis models. 

The scale height of subdwarfs is rarely discussed. We 
used the two samples here to estimate the scale height at 
l±0.5kpc from the peak in their distance distributions near 
2± lkpc (Figure[5]). However, to do so we must assume that 
the each subdwarf plus companion system (independent of 
system parameters) is drawn from the same parent distance 
distribution, and we pick each of these with the same fre- 
quency. If the scale height is ~ lkpc, it is therefore most 
consistent with the Galactic thick disk scale height (e.g. 
0.75 ± 0.07kpc. IdeJong et allboiol ). If the subdwarf pop- 
ulation was associated with the th in disk, a smalle r scale 
height of 0.3 kpc would be expected (|Juric et a l. 2008), while 
a rise towards 25 kpc would have indicated a halo population 
|de Jong et al.ll2010l ) . More accurate modelling of individual 
subdwarfs together with a larger volume limited sample is 
required to study the distribution and reliably quantify the 
scale height of the subdwarf population. Our methods are 
well suited to offer such large samples as ongoing and near- 
future surveys cover an increasing part of the sky. 



9 CONCLUSION 

We have developed a method to select hot subdwarfs stars 
with mid-M to early-F-type near main-sequence compan- 
ions using a combination of ultraviolet, optical and infrared 
photometry. This selects a complementary sample to those 
found from radial velocity surveys, which typically limit 
themselves to objects with no obvious evidence for a com- 
panion in the optical range. We applied this method to two 
samples, one selected from a match between GALEX, CMC 
and 2MASS (covering a large area), and the other using 
GALEX, SDSS and UKIDSS (probing deeper in the K 3 - 
band and therefore further away). We also use the SDSS for 
fitting in the C2MS sample. 

A significant number of subdwarfs with F to K-type 
companions were found in both samples. The distributions 
are consistent with the systems being produced, at least 
in a significant part , by the very efficient RLOF chan- 
nel (lHan et al. 20031 ) . However, nei t her th e predictions of 
lHan et alj (120031 ) or IClausen et~ai1 (|2012l ) match the ob- 



served distribution completely. We find that M-type com- 
panions are far less prevalent than K-type systems. 

It is clear that, at least for a large fraction of the subd- 
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warf population, prior binary evolution plays an important 
role. This group has largely gone unstudied previously. With 
future surveys such as the Southern SkyMapper project and 
VISTA, the same procedure as carried out here can be ap- 
plied to a large field in the southern sky. This would find 
many more subdwarfs with early type companions and allow 
for a thorough test of our understanding of the prior binary 
evolutionary pathways required to form the large subdwarf 
populations we see. Similarly, the Wide-field Infrared Sur- 
vey Explorer could be an excellent addition to this search, 
allowing us to probe for fainter companions and covering the 
whole sky. 
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Table 9. Example of: Full list of objects from the C2M sample with magnitudes inside the cuts described in TableTT] Online-only Table. 
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Table 10. Subdwarf and companion effective temperatures, and distance estimates for the C2MS sample of subdwarf plus companion star systems when fitted using the method 
described in Section[6] A MST G ff = OK corresponds to the best fit being a single subdwarf, without the addition of a companion. The "{" notation is described in Section[7] and does 
not simply repre sent uncertainties. In all cases, a more realistic error on the subdwarf temperatures is a few thousand Kelvin (see Section lY.^H . The "E(B-V)" column is the reddening 
according to the lSchlegel et ah I Jl998ft map. The "Q" (Quality) column values correspond to; l:Good fit, 2:Average fit, 3:Poor fit, 4:WD/WD+MS/CV and 5:Quasar/Galaxy. Values of 
three and above are excluded from the histograms shown in FigurelBl l9l and fTOl The classifications in this catagory between values of 1, 2 and 3 are purely qualitative. The "SDSS Spec" 
column cor responds to the visua l classific ation of the SD SS spectrum. The "K no wn Comp" highlights ob jects wh ich were known previously to be com p osite subdwarf plu s c ompanion 



syste ms (l:lFerguson et al.l 1984 2: iKilkennv et al.lll988l. 3: lAllard et aTJll994l, 4: lAznar Cuadrado fc Jeffery||200ll. 5: lAznar Cuadrado fc Jeffery||2002l, 6: iLisker et al]|2005l. 7: ICstensen 
l2006l. 8: IStroeeret*^l*ll2*00*7 ). The final column shows objects also included in TablelBI where "pWD" stands for possible white dwarf, according to the classification of lGirven^t**al 
(2011). 
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07:10:29.29 


+29:38:52.3 


21/22 


6 75U-9S 

u ' ' 16.50 


1 5/?i 
x - u 11.5 


25/25 

^124 


6 75(1-22 

u - 1 J 16.50 


1 5((2 
X - J ll.3 


0.074 




1 






0735+2012 




07:35:46.24 


+20:12:35.6 


21/22 


5 5015-15 

J,,JU 15.25 


2 l/3 f 


23/2f 

^-J 121 


5 25/5-52 
"•""IB. 00 


1 9(?-i 


0.041 




2 






0737+2642 




07:37:12.24 


+26:42:25.3 


25/26 

J 124 


5 50/5' 15 


1 6(4-f 


27/21 
* ' 126 


5 2515-52 


1 5(H 

± . J 11.4 


0.039 


WD 


1 


SD 




0754+1822 




07:54:04.24 


+18:22:40.4 


22( 23 


7 25( 7 -52 


3 4/3- 7 

J -^13.1 


25/26 

^ u 124 


1 - ,J1J l7.25 


3 5(2-2 


0.045 




1 






0755+2128 




07:55:49.51 


+21:28:18.0 


17/18 
x ' 116 


7 25/!'5R 


2 3/3 g 

J 12.3 


20/21 


7 75(1-22 

' • ' '-'1-7.50 


2 5(2-5 


0.065 




1 






0804+2250 




08:04:20.93 


+22:50:18.0 


37(42 

Jl 136 


5 75(1-2°, 

° ' J 15.50 


1 5fH 

1,J ll. 5 


31122 


4 50(1-15 


1 o(i-2 

1,U 10. 9 


0.048 




3 






0805—0741 




08:05:16.32 


—07:41:50.6 


^ a 128 

14(45 
"113 


7 oo/Z'25 


2 4/2.8 

12.2 


28(2? 


7 oo(I-25 

1 - uu l-6.75 


2 I/? ? 


0.117 




1 






0812+1911 




08:12:56.86 


+19:11:57.9 


7 25( 7 -52 


^•°14.8 


15/16 


7 oo(I-2 5 


4 4(4-2 


0.035 


CV 


4 


CV 


CV 


0814+2019 




08:14:06.84 


+20:19:01.0 


21/22 


5 50/5 15 
J " JU 15. 25 


2 0/3-2 
^■ U 12.0 


23/24 

^ u 121 


5 50(5-15 

a - uu l5.25 


2 1(2-2 
^■■ L 11.8 


0.042 




1 


SD 




0815+4740 


PC 0812+478 


08:15:48.88 


+47:40:40.4 


40(42 


5 OOI5 25 

UiW l4.75 


1 4/1-1 

■""^l 1.3 


40(42 
w 139 


4 75(5-22 

' " J 14.50 


1 2(}2 


0.067 


WD 


2 






0818—0701 




08:18:06.86 


—07:01:23.9 


22/2? 
^121 


7 75/I-9S 

' - ' J 17.50 


3 4-fH 
J -" 3.1 


14/15 
"113 


5 75(2-22 

u - ' u 15.50 


2 4(2-2 

3. 4 


0.097 




1 






0820+1739 




08:20:03.34 


+17:39:14.0 


20(?i 
Z,U 119 


6 75/I go 

u ' J 16.50 


2 9/2-9 

*' n 11.8 


20i?i 


6 50(2-15 

U " JU 16.25 


2 6(2-2 
^■"11.6 


0.033 




1 


SD 




0824+3028 


PG 0821+306 


08:24:34.03 


+30:28:54.6 


2i{i 


5-25{ 5 ; 00 


i-4(?:l 


21(1§ 


5.00{|;2| 


1-3{?| 


0.044 




1 


SD 




0825+2006 




08:25:07.22 


+20:06:36.5 


24{|I 


6.00{|| 




27(i 


6-oo{|:i 


1-7(1:1 


0.037 




1 






0825+1202 




08:25:44.73 


+12:02:45.2 


22{J? 


8.25{ 8 : o 


2-5(2!3 


26{^ 


8.75(1:88 


2-7(2.8 


0.044 




1 






0825+1307 




08:25:56.86 


+13:07:54.3 


24{|I 


5.00{!J 


i-o(ii 


27{|| 


5.00{| : 2| 


i-i(li 


0.034 




2 






0829+2246 




08:29:02.64 


+22:46:37.6 


26{|I 


6-00{|;2| 


2-7(li 


28{i 


6.oo{|;| 


2.7{|| 


0.036 




1 


SD 




0833-0006 




08:33:37.88 


-00:06:21.4 


29{|g 


7 25( 7 - 50 
' - zo l7.00 


3-1(1:1 


29{i 


7.00{ 6 7 ;i 


2.5(2:3 


0.041 




2 






0844+3102 


PG 0841+312 


08:44:08.18 


+31:02:09.3 


22{|f 


4.75{t;go 


i.o{J:5 


22{|3 


4.75(i88 




0.049 




1 






0849+1337 




08:49:51.40 


+13:37:00.4 






2-4{| : ! 


21 /22 


6.50{| : || 


2.0(|| 


0.040 




2 






0907+2739 




09:07:34.26 


+27:39:03.4 


2i(i 


5-50{l;i 


2-5{|;g 


21(1 


5-50{l; 7 | 


2.3(2:3 


0.026 


WD 


3 






0923+0652 




09:23:58.62 


+06:52:18.3 


29{i 


6-75{J:88 


2.i(l:§ 


29(lg 


D . zo 16.00 


i-5(?:l 


0.054 




1 




pWD 


0924+2035 


PG 0921+208 


09:24:05.20 


+20:35:46.8 


19(?g 


4-25{f:|g 


1-6(1:1 


19{?g 


3.00{g;8g 


i.edl 


0.041 




3 






0929+0603 




09:29:20.48 


+06:03:47.1 


2H18 


5-75{l;g8 


i-6(?:I 


29(ig 


5.50{g;H 


i-2(?:8 


0.052 




2 






0935+1621 


PG 0932+166 


09:35:41.37 


+16:21:11.0 


30{| 


5.00{f;| 


i-o(S:g 


29{| 


4.75{|;88 




0.033 




1 






0937+0813 


PG 0935+084 


09:37:40.95 


+08:13:20.5 


23(22 


6-00{«;?l 


2.0{2;| 


23{ 2 2 


5-75(1:88 


mii 


0.042 


sdB 


1 


SD 




0941+0657 


PG 0939+072 


09:41:59.35 


+06:57:17.2 


2i(i 


6-25{i:gg 


i-7(?:I 


21(i 


6.00{| : 2| 




0.040 


WD 


1 






0958+2236 




09:58:15.97 


+22:36:04.2 


33{|f 




1.8{2;0 


34(I| 


6.75{ 6 7 :88 


i-8(?:? 


0.033 




1 






1003+3716 


PG 1000+375 


10:03:19.69 


+37:16:35.1 


30(31 


5-50{|;5| 


i-2(i:i 


29{i 


c w r5.50 
°. z °l-5.00 


1.1(1:1 


0.016 


WD 


1 


SD 




1005+4317 




10:05:05.07 


+43:17:36.5 


29(i 


6-75{S;gg 


i-9(?:I 


29{i 


6.75{ 6 7 :88 


i. 9 (f:l 


0.012 




1 






1015-0308 


SW Sex 


10:15:09.39 


-03:08:32.3 


18(1? 


7-00{^| 


2-6(1:1 


mil 


7.00{ 6 7 : 2| 


2.6{2 : 6 


0.033 


NL 


1 




CV 



to 

Table 10 - continued >t^ 



Name 


Identifier 


R.A. 






Dec 








1018+0721 




10: 


:18: 


:01. 


.55 


+07: 


21: 


24. 


.4 


1018+0953 




10 


:18: 


33 


.15 


+09: 


:53: 


36. 





1027+2409 


PG 1025+244 


10 


:27 


:51 


.19 


+24: 


09: 


17. 





1049+1842 


PG 1046+189 


10 


:49 


:33 


.53 


+ 18: 


12: 


11. 


.5 


1100-2113 


EC 10583-2057 


11 


:00 


:46. 


.69 


-21: 


13: 


12. 


.3 


1102+2616 




11 


:02 


:11 


.09 


+26: 


16: 


16. 


.3 


1113+0413 


PG 1110+045 


11 


:13 


:17 


.31 


+04: 


13: 


11. 


.7 


1131+0932 


PG 1128+098 


11 


:31 


:14 


.37 


+09: 


32: 


20. 


A 


1149+2231 


PG 1146+228 


11 


:19 


:00 


.50 


+22: 


31: 


05. 


.9 


1203+0909 


PG 1200+094 


12 


:03 


:19. 


.46 


+09: 


09: 


51. 


6 


1212+4240 


PG 1210+429 


12 


:12 


:38 


.56 


+42: 


10: 


02. 


1 


1233+0834 




12 


:33: 


:09 


.68 


+08: 


31: 


31. 


1 


1316+4359 


PG 1314+442 


13: 


:16 


:33 


.00 


+43: 


59: 


01. 


.9 


1325+1212 


PG 1323+125 


13: 


:25 


:57. 


.21 


+ 12: 


12: 


20. 


6 


1326+0357 


PG 1323+042 


13: 


:26 


:19. 


.95 


+03: 


:57: 


51. 


.3 


1402+3215 




14 


:02 


:32 


.86 


+32: 


15: 


21. 


5 


1404+2450 


PG 1402+251 


14 


:04 


:29 


.98 


+24: 


50: 


20. 


6 


1407+3103 




14: 


:07 


:47. 


.63 


+31: 


03: 


:18. 


.3 


1421+0753 


KNBoo 


14: 


:21 


:38 


.21 


+07: 


:53: 


20. 


.9 


1502-0245 


PG 1459-026 


15 


:02 


:12 


.13 


-02: 


15: 


56. 


.7 


1517+0310 


PG 1514+034 


15 


:17 


:14 


.30 


+03: 


10: 


27. 


6 


1518+2019 


PG 1516+205 


15: 


:18 


:38 


.81 


+20: 


19: 


17. 


.0 


1524+0134 




15: 


:24 


:03 


.04 


+01: 


31: 


21. 


.3 


1528+1300 




15 


:28 


:33 


.90 


+ 13: 


00: 


:57. 


2 


1530+1204 




15 


:3Q 


:05 


.00 


+ 12: 


01: 


02. 





1542+0056 




15: 


:42 


:18. 


.31 


+00: 


56: 


12. 


6 


1602+0725 


PG 1559+076 


16: 


:02 


:08. 


.96 


+07: 


:25: 


10. 


.8 


1603+0954 




16 


:03 


:24 


.51 


+09: 


:54: 


.12. 


,9 


1610+3450 




16 


:10 


:40 


.72 


+34: 


50: 


.11. 


1 


1618+2141 




16: 


:18 


:06 


.46 


+21: 


11: 


25. 





1619+1453 


PG 1617+150 


16: 


:19 


:49. 


.35 


+ 14: 


:53: 


09. 


.5 


1629+0026 




16 


:29 


:06 


.76 


+00: 


:26: 


19. 


.9 


1640+3842 


PNG061.9+41.3 


16 


:4Q 


:18 


.20 


+38: 


12: 


20. 


6 


1709+4054 


PG 1708+409 


17: 


:09 


:59 


.23 


+40: 


51: 


19. 


.5 


1710+2238 




17: 


:10 


:36 


.45 


+22: 


:38: 


07. 


1 


1734+3213 




17 


:34 


:49 


.30 


+32: 


13: 


13. 


.5 


1822+4320 




18 


:22 


:42 


.87 


+43: 


20: 


37. 


.4 


1 8344-49*37 




18: 


:34 


:14 


16 


4-49' 


37 


27 




2020+0704 




20 


:20 


:27. 


.24 


+07: 


01: 


14. 


.5 


2023+1230 




20 


:23 


:14 


.11 


+ 12: 


30: 


56. 


.9 


2047-0542 




20 


:47 


:42 


.36 


-05: 


12: 


32. 





2052-0457 




20 


:52 


:26 


.19 


-04: 


:57: 


16. 


.0 


2056+0425 




20 


:56 


:19. 


.33 


+01: 


:25: 


23. 


6 


2105+1635 




21 


:05 


:15 


.38 


+ 16: 


35: 


:18. 


3 


2117-0015 




21 


:17 


:15 


.90 


-00: 


15: 


17. 


.7 


2117-0006 




21 


:17 


:42. 


.22 


-00: 


06: 


19. 


.9 



No Correction 
sdB T eff MS T ea d 
(1000 K) (1000 K) (kpc) 



Reddening Corrected 



sdB T eff 


MS T eff 


d 








SDSS 


(1000 K) 


(1000 K) 


(kpc) 


E(B-V) 


SIMBAD 


Q 


Spec 


30{i 


5-75{f;gg 


i-6{i:I 


0.027 




3 


SD 


29{i 


4.50{|;| 


0.8{i;| 


0.037 


WD 


1 




26{i 


6-25{ 6 ; 00 


2-i{?:S 


0.017 




1 


SD 


20{| 


5.50{|; 7 1 


2.0{?;0 


0.033 




1 






5.75{6;» 


i-4{J:l 


0.053 




1 




2Hi 


6.00{f ;7 | 


i.7{? : ? 


0.019 




1 


SD 


29(1 


4.50{f; 7 l 


o.8{5 : i 


0.051 




1 




40{te 


5.75{g;g8 




0.039 




2 


SD 


24{| 


5-25{|;ES 


1-4{H 


0.022 




1 


SD 


28{ 2 6 


5-75{l;gg 


i-5{J:8 


0.020 




1 




24{|I 


5.75{6 : « 




0.015 




1 


SD 


29{1 


5-75{|;|g 


i-6{?:l 


0.019 




2 


SD 


28{|? 


5.00{f;| 


i-e{i:§ 


0.021 




1 




27{i 


5.75{g;ge 


2-i{?:l 


0.034 




1 




26{| 5 7 


5.00{|;?1 




0.025 


sdO 


2 


SD 


22{|? 


6.00{|;?I 


1-7{J;| 


0.015 




1 


SD 


27{I 


6.00{|;| 




0.017 




1 




mil 


4.75{f;g» 


i-4{J:l 


0.011 




1 




29(1 


4.50{|;5| 


0.9{i;| 


0.028 


sdB 


1 


SD 


24{i 


6-00{|;?I 


i-s{i:I 


0.124 




1 


SD 


39{i 


5.75{|;°° 


0.9{J : o 


0.039 


WD 


1 




26{| 5 7 


6-00{6;?| 


i-5{};I 


0.051 




1 




28{|? 


6-00{«;?I 


i-i{J:J 


0.065 




1 




40{|g 


5-75{|;g8 




0.040 




3 




14{?J 


9.25{|;| 


4.8{S : l 


0.038 




2 




2Ui 


6-25{ 6 ;o 


i-3{f:§ 


0.098 




1 




37{| 


5-25{ 5 ;q 


i-o{J:§ 


0.047 




1 




32{2g 


5-75{«;g° 


2.0{?; 7 


0.059 




3 


SD 


31{p 


7-50{H| 


2-o{!:5 


0.018 




1 




24{|f 


6-00{I;i 


i-5{J:3 


0.069 




1 




29{i 


5-75{|;| 




0.051 




1 




32{i 


7-25{ 7 ;o 


2-4-Cl;| 


0.095 




1 




16{1 5 7 


3.00{|2S 


1.8{i;f 


0.012 


PN 


3 




28{|? 


5-50{g;5g 


1.8ft? 


0.029 


WD 


1 




34{H 


7.50{ 7 ; 7 | 


i-4{l:t 


0.058 




1 




25{H 


6.00{|;?1 




0.054 




1 




40{|g 


o.oo{% 00 00 


0.6{J° 6 7 


0.050 




1 




23{| 


6.00{6;2| 


2.3{l;f 


0.055 




3 




38{ 37 


4.50{|;7S 


0.6{°;| 


0.143 




1 




35{|I 


4.75{t;« 


0.8{0;B 


0.123 




2 




24{|i 


5-75{6;» 


i-4{}:I 


0.050 




1 




37{i 


6.25{g : gg 


i-3{i:§ 


0.104 




1 




21 /22 
zi t20 


6-00{I; 2 7 I 


2.2{|? 


0.092 




1 




30{i 


4.75{| : oo 


0.8{0;1 


0.075 




3 




14{?§ 


12.75{13;00 


5-4{l| 


0.055 




1 




22{|? 


6.00{6;2| 




0.074 




1 


SD 



re 



29{|g 
28{|? 



"123 

mil 

30{3i 
22{2f 
29{| 
38{ 37 
23{i 
27{I 
23{i 
30{i 
28{i 

mil 

24{i 
22{i 
27{I 

mil 

27{I 
24{il 

{40 
M!) 



-..7", 

(i.oo 



40{1 
24{i 
27{I 
40{|§ 

IHli 
29{| 

32{3o 
29{i 
29{i 
22{2f 

29{i 
29{1 

mil 

31{i 
23{i 
3Hi 
22{| 

23{ 2 2 
2U 22 

22{i 
29{i 
22{2f 

21{i 

mil 



6.00{ 625 
5.75{ 
6.Q0{§; 7 ° 

5- 75{I:i» 
6.25{t : |j 

6- 25{ 6 ; o 
4.75{| : oo 
5-75{g;gg 
5.25{| : Sg 
5.75{|;|g 
5-75{|:ig 
6.00{I;f| 

5.25{|:S 

5- 75{| : |g 
5.00{|; 2 1 

6- 25{|;g§ 
6.25{|:§8 
5.25{|| 

6-25{|;g° 
6.00{l : f| 
6.00{|;| 

6.oo{I : f| 

8-00{f;18 

6- 50{l; 7 l 
5.25{§ : §8 
5-75{I:ig 

7- 50{ 7 ; 7 | 
6.00{|;|1 
6.00{|;|| 
7 25l 7 - 50 
3-00{|;gg 

5- 50{p| 
7 25l 7 - 50 
6.00{|;2| 
3-00{|;gg 
6.25{| : Sg 

5.75{| : g8 
6.00{| : 2| 

6.00{I;f| 

6- 00{|:f| 
6-50{|; 7 l 
6.00{| : 2| 



1.9{? : ° 
1.9{f : | 

2-i{?:J 
2-o{?:l 

0.9{ ;| 
1.2{i;f 

i-4{y 

1- 5{l:I 
i.9{? :7 

2- i{?:| 
i-5{J:S 
i.9{f :7 

l-8{?;g 

i-5{J:S 

m\% 

i-5{l:I 

1- 3{1;1 
4.8{|;| 

0.9{ ;i 

2- o{|;g 

2.0{1;3 

i-3{l:I 

2.5{1;° 

1- 4{}:I 
o.5{S!:I 

2- 7{l;§ 

i.7{? : | 
i-6{}:! 
i-2{l:I 

3- 0{|; 7 

1- 8{?:| 

3.5{ 2 ; 3 

2- 0{?;| 



2,7 



1,2,7 SD 



1,7 SD 



SD 



pWD 
pWD 



® 



Table 10 - continued 



No Correction 



Reddening Corrected 



Name 


Identifier 


R.A. 


Dec 


sdB T cff 
(1000 K) 


MS T eff 
(1000 K) 


d 

(kpc) 


sdB T cff 
(1000 K) 


MS T eff 
(1000 K) 


d 

(kpc) 


E(B-V) 


SIMBAD Q 


SDSS 
Spec 


Known 
Comp 


Tabled 


2129+0045 




21:29:06.05 


+00:45:09.6 


25{|I 


5.50{i; 2 ! 


2-o{?;l 


29{|| 


4.75{i : §° 


l.2{{ 


.9 
.2 


0.044 


WD 


1 


SD 










91 -9Q-9Q 1 1 
z_L.zy.zy.ir 


4-1 O^Q-OQ Q 


91 / 2 3 


°- ou l5.25 


1 5/ 24 


zo \24 


t; c n /5.75 
°- ou i5.25 


1.5{} 


.(. 
.4 






1 








2135+2026 




21:35:51.03 


+20:26:45.3 


21{i 




i-8{?:I 


22{|3 


5.00{f;i 


1A{\ 


.5 
3 


0.124 




1 








2138+0442 


PC 2135+045 


21:38:00.77 


+04:42:11.5 


25{^ 


5.00{|;?§ 


1.2ft? 


27{|I 


4.oo{|:i 


1.0{\ 


1 

.0 


0.056 




1 




3,4,5 


SD 


2143+1244 




21:43:54.65 


+12:44:58.3 


30{3i 


6-25{|;gg 


2-9{l;l 


29(1 


5-75{g : S8 


2-1(1 


1 
.0 


0.095 


CV 


4 


CV 




CV 


2147-0837 




21:47:08.07 


-08:37:47.5 


29{|g 


6-50{|il 


2-3{| : ; 


29{|? 


6 25l 6 - 50 
D - zo l-6.00 


1.9{? 


.2 

s 


0.044 




1 


SD 






2223+3850 




22:23:26.78 


+38:50:16.7 


32{I? 


7.75{B ; oo 


2-9{Ii 


28{|9 


8.00{?;i 


2-8(1 


1) 
.2 


0.108 




1 








2346+3657 




23:46:21.39 


+36:57:27.6 


29{|g 


7 25l 7 ' 50 


2.5{§;2 


22{|f 


6.75{|;gg 


2-0{f 


2 

.K 


0.142 




2 








2346+0344 




23:46:55.71 


+03:44:29.4 


30{3i 




i-ni:! 


29{30 


6-50{g : ™ 


1.5{2 


2 

.4 


0.054 




1 
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Table 11. Subdwarf and companion effective temperatures, and distance estimates for the SU sample of subdwarf plus companion star systems when fitted using the method described 
in Section[6] following the same format as Table lTOl A MS T c fj = K corresponds to the best fit being a single subdwarf, without the addition of a companion. The "{" notation is 
described in Section]?] and does not simply represent uncertainties. In all cases, a more realistic error on the subdwarf temperatures is a few thousand Kelvin (see Section[7TH]l . The final 
column shows objects also included in Table[6] 



No Correction 



Reddening Corrected 











sdB T a 

bUD ± e ff 


MS T a 




sdB T a 


MS T a 


( 1 








SDSS 


1\.11U w 11 




Ntirnc 


Identifier 


R.A. 


Dec 


C1000 TO 


("1000 JO 


Ckncl 


nooo TO 


nooo 10 




E(B-V) 


SIMBAD 


O 


Spec 


Comp 


Table[6] 


0018+0101 


HE 0016+0044 


00:18:43.52 


+01:01:23.6 


24{i 


4.25{|;gB 


i-2{i:I 


28{| 7 


4.75{i;0° 




0.029 


sdB 


1 






SD 


0032+0739 


PB 6015 


00:32:21.87 


+07:39:34.4 


21{! 


5-25{g-gg 


4.4{f;l 


21{i 


5.00{|?| 


4.1{|f 


0.040 


Comp 


3 






WDMS 


0051+0059 




00:51:49.65 


+00:59:50.7 


27{|I 


5.00{|;| 


4.5{1;| 


28{| 7 


5.00{|;| 


4-6{f:l 


0.028 


WD 


5 








0054+1508 




00:54:11.11 


+15:08:19.3 


29{i 


7.25{ 7 ; 00 


2-8{|| 


29{30 


7.25{ 7 ; 00 


2-7{|:I 


0.059 




3 








0141+0614 


HS 0139+0559 


01:41:39.92 


+06:14:37.6 


14{i| 


7.25{ 7 ; g 


3.7{f; 7 


14{H 


7-00{|;?I 




0.048 


NL 


2 






CV 


0205+0712 


PB 6645 


02:05:59.87 


+07:12:39.9 


25{i 


5-75{|;gg 


7-3{|;| 


Mil 


5-75{f;°° 


7-3{I;| 


0.056 




1 








0300-0023 


WD 0257-005 


03:00:24.57 


-00:23:42.1 


38{|! 


5.00{tf B 


3-2{2i8 


39{|g 


4.50{£™ 


2-6{|? 


0.120 


Comp 


1 






WDMS 


0316+0042 


PG 0313+005 


03:16:20.13 


+00:42:22.8 


Mil 


6-00{|;?l 


2-o{?;I 


27{|I 


6.00{6;?5 


2-o{ 2 : § 


0.087 


WD 


1 


SD 






0737+2642 




07:37:12.27 


+26:42:24.7 


MH 


5-50{lil 


i-«Kl:§ 


26{i 


5.00{| : 2 1 


i-4{i:i 


0.039 


WD 


1 


SD 






0744+2103 




07:44:41.80 


+21:03:52.6 




4.75{l;go 


2.8{§;° 


39{|g 


4-00{f; 2 7 | 


2-3{l;I 


0.053 




1 








0752+2535 




07:52:39.82 


+25:35:50.0 


mil 


4.25{igg 


4-8{I;| 


"{15 


18.00{i 7 ;°o 

e-MHI 


36.9{i;l 


0.061 


WD 


2 








0755+2128 




07:55:49.49 


+21:28:18.5 


21{ il 




2-o{? : | 


"{II 


i-7{?:| 


0.065 




1 








0758+2818 




07:58:13.60 


+28:18:16.0 


23(22 


6-50{|il 


3.6{|;1 


18{« 


6 25{|;gg 


3.5{|;? 


0.039 




1 


SD 






0809+1924 




08:09:21.89 


+19:24:00.1 




°- 00{ |°lo 


2.3{| 6 2 6 


40{|g 


0-00{§5°8o 


2-0{? 1 9 4 


0.038 




3 


WD 






0813+2813 




08:13:52.02 


+28:13:17.3 


Mil 


6-25{e!oo 


6A{ ii 


Mil 


6.00{1;?1 


5.8(2:S 


0.032 


cv 


4 


CV 




CV 


0814+2019 




08:14:06.80 


+20:19:01.7 


M% 


6-25{||g 


3.5{ 27 


Mil 


5-75{g;°° 


2.9{|} 


0.042 




1 


SD 






0814+2811 




08:14:53.92 


+28:11:22.5 


22{|? 


6.00{|| 


3-7{t:i 


Mil 


5.75{|;| 


3.2{ 2 ;g 


0.030 




1 


SD 




pWD 


0829+2246 




08:29:02.62 


+22:46:36.8 


21{! 


5-25{ 5 ;oo 


i-9{?:I 


Mil 


^■25{5!oo 


1.9{?;° 


0.036 




1 


SD 






0833-0006 




08:33:37.87 


-00:06:21.4 


29{1 


^Hii 


2-4{l:I 


Mii 


6.75{^° 


2-3{| : f 


0.041 




2 








0843-0048 




08:43:51.07 


-00:48:24.6 


Mil 
24{|1 




3.o{| : § 


19{fg 


6.00{|;?I 


2.7{|;2 


0.033 




1 








0846+0142 




08:46:28.66 


+01:42:17.0 


4-50{i^l 


4.6{£§ 


Mfr 


4.00{|; 2 1 


4.5{f; 7 


0.042 




1 








0854+0853 


PNA66 31 


08:54:13.16 


+08:53:52.9 


40{|B 


3.00{|| 


i-2{i:i 


Mf 9 




l.l{^i 8 


0.065 


PN 


3 






pWD 


0856+0518 




08:56:33.17 


+05:18:39.6 


29{1 


3-25{ 3 ; 00 


2-6{f:l 


4 ° { 1 




3-l{|| 


0.050 




1 








0859+0759 




08:59:26.08 


+07:59:13.3 


22{|f 


5.75{|;gg 


2.9{|| 


23(22 


5.50{|| 


^■^{2.3 


0.080 




1 








0902+0734 




09:02:25.06 


+07:34:04.0 


21{i 


6-50{Iil 


3.5{ 33 


22{|f 


6-25{e!oo 


3-i{i:I 


0.072 




1 


SD 






0906+0251 




09:06:40.00 


+02:51:46.4 


39{|i 




5-7{|;i 


Mfr 


6.50{|| 


5-2{|;6 


0.034 




1 


SD 






0906+0437 




09:06:00.86 


+04:37:45.1 


21{! 
34{|1 


5.75{g : gg 


2-4{|;| 


Mil 


6-25{ B ; o 


3.0{|;| 


0.037 




1 


SD 






0920+1057 




09:20:48.04 


+10:57:34.5 


4.75{t;0[J 


3.3{|f 


40{« 


4.75{t;g[J 


3.3{ 3 ; 2 


0.036 


Comp 


1 






WDMS 


0920+3356 


BKLyn 


09:20:11.21 


+33:56:42.4 


Mil 


675{||g 


2.3{?; 7 


20{?J 


6-50{|;ll 


2-l{|;| 


0.017 


NL 


3 


WD 




CV 


0925-0140 




09:25:35.00 


-01:40:46.8 


mil 


5-50{li| 


9.4{J 4 / 


Mil 


5-50{Ii| 


9-3{|;| 


0.031 




2 








0929+0603 




09:29:20.43 


+06:03:46.2 


29{i 


6.00{f;?l 




29{| 


5-50{§ii 


1.2{l;f 


0.052 




1 


SD 






0937+0813 


PG 0935+084 


09:37:40.93 


+08:13:20.9 


21{|g 


5-75{|;g8 


i-7{?:l 


23{22 


5-75{g;°° 


l-7{i| 


0.042 


sdB 


1 


SD 






0939+3038 




09:39:14.38 


+30:38:17.3 


39{f 3 


6.25{f; 7 ° 


5.0{I;f 


37{H 


6.00{|| 


4-4{|:| 


0.017 




1 


SD 






0941+0657 


PG 0939+072 


09:41:59.32 


+06:57:17.2 


29{li 


6-50{|il 


i-n?:g 


30{3i 


6-25{ eoo 


i-5{l:I 


0.040 




2 








0951+0347 




09:51:01.29 


+03:47:57.0 


23{i 


4.00{|;?| 


1.9{?;° 


28{| 7 


4-00{|; 2 l 


2.0{ 2 ;0 


0.039 


WD 


1 


SD 






0959+0330 


PG 0957+037 


09:59:52.01 


+03:30:32.8 


31{i 


3-00{|;25 
5.00QM 




40{fg 


; ;V: 




0.025 




2 






pWD 


1006+0032 


PG 1004+008 


10:06:45.75 


+00:32:04.5 


Mil 


3.3{|;f 


29{i 


4.25{tgg 


i-8{?:J 


0.036 




1 






pWD 


1011-0212 




10:11:36.23 


-02:12:14.6 




4.75{i;go 


4.0{f; 7 


28{| 7 


4.75{t;00 


4.2{f : | 


0.040 




1 








1012+0044 




10:12:18.95 


+00:44:13.4 


Mil 


4.75{f : °o 


5.0{|;| 


29{| 


3-75{|;°o 


2.9{f;l 


0.033 


WD 


1 


SD 






1015-0308 


SW Sex 


10:15:09.38 


-03:08:32.8 






2-o{? : | 


23(22 


6-50{|;5I 


mii 


0.033 


NL 


2 






CV 


1016+0443 




10:16:42.94 


+04:43:17.7 


29{i 


4.50{*™ 


4-8{I;l 


29{| 


4.00{|;| 


4-3{Ii 


0.024 


Comp 


i 


WD 






1018+0953 




10:18:33.11 


+09:53:36.1 


35{|1 


5-50{|il 


i-3{i:l 


38(37 


5-25{ 5 ; 00 


1.2{};f 


0.037 


WD 


1 


WD 






1034+0327 


HS 1031+0343 


10:34:30.16 


+03:27:36.4 


18{« 


4-25{t§8 


4.7{|;° 


20{« 


4.50{« 2 g° 


5.0{|y 


0.038 


WD 


3 
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Table 11 - continued 



Name 


Identifier 




R.A. 






Dec 








sdB T off 
(1000 K) 


MS T eff 
(1000 K) 


d 

(kpc) 


sdB T eff 
(1000 K) 


MS T eff 
(1000 K) 


d 

(kpc) 


E(B-V) 


SIMBAD 


Q 


SDSS 
Spec 


Known 
Comp 


Table[6] 


1040+0217 






10: 


10: 


32 


.74 


+02: 


:17 


:29 


.7 


30{1 


4.25{|;g8 


2.4(3:1 


39{|g 


4.50{| : 7| 


3.0{| 


1 


0. 


.035 




1 








1055+0930 






10: 


.55: 


:25. 


.88 


+09 


:30 


:56 


.3 


28{|9 


6-50{|;i 


mil 


27{|I 


B.oo{g : « 


5.5{| 


!) 


0. 


.030 




1 


SD 






1057-0230 






10: 


:57: 


59. 


29 


-02: 


:30: 


:02 


.1 


27{|1 


5.50{|;| 




27{|I 


5-25{f : §H 


5.1{| 


4 
.7 





.042 




2 








1100+0346 






11: 


00: 


53. 


55 


+03: 


:46: 


:22 


.8 


34{|I 


3-75{ 3 ; 5 o 


2.7{I;1 


40{|g 


3-75{|;gg 


2.9{| 



!' 





.044 


WD 


1 


SD 




pWD 


1113+0413 


PG 11104 


-045 


11: 


13: 


17. 


32 


+04: 


:13: 


:14. 


.5 


30{3i 


4.75{f : o° 


o.9{J : g 


34{|| 


5.00{f : 2 l 


1.0{1 


1 
!' 


0. 


.051 




1 




2,7 




1116+0755 






11: 


16: 


16. 


37 


+07: 


:55: 


:32. 


.5 


28{|? 


5.00{t;?| 


2-3{?:i 


29{i 


3-75{f:°8 


1.2{i 


!) 


0. 


.042 




1 






pWD 


1131+0932 


PG 1128+098 


11: 


31: 


14. 


32 


+09: 


:32: 


:19. 


.0 




5-75{f;°° 




39{fg 


5-50{l: 2 7 I 


1.1{1 


; 
(1 





.039 




1 








1134+0153 






11: 


3d: 


18. 


00 


+01: 


:53: 


:22 


.1 


38{1? 


6.50{g;?| 


mil 


24{|I 


6 ,oo{l: 2 l 


5.8{| 


r . 





.031 


WD 


1 


SD 






1135+0731 






11: 


35: 


36. 


86 


+07: 


:31: 


:28. 


.3 


29{i 




6.4{|;f 


29{|g 


6.00{1; 2 1 


5.5{T 

1-4. 


3 

:) 


0. 


.041 




1 


SD 




pWD 


1203+0909 


PG 12004 


-094 


12: 


03: 


19. 


38 


+09: 


:09: 


:51 


.6 


27{|| 


5.75{i : gg 


i.5{i:I 


28{|? 


5-75(1:88 


1.5{1 


'> 


0. 


.020 




1 








1215+1351 






12: 


15: 


23. 


73 


+ 13: 


:51 


:02 


.3 


mil 


4.50{*;™ 


3-i{|:l 


21{ 1 


4.oo{| : || 


2.9{| 


8 





.032 




2 


SD 




pWD 


1228+1040 


WD 1226- 


f 110 


12: 


28: 


59. 


93 


+ 10: 


: 10: 


33 


.0 


21{! 


3-00{l;gg 


2.2{|:I 


23(22 


3-oo{l : gg 


2.3{| 


4 





.028 


WD 


4 


WD 




pWD 


1233+0834 






12: 


33: 


09. 


62 


+08: 


:34: 


:3 d 


.5 


30{| 


6.00{|;2| 


i.fl{?:? 


30{3i 


5-75{| : 9g 


1.7{1 : 


n 

.5 


0. 


.019 




1 


SD 






1235+1029 






12: 


35: 


13. 


03 


+ 10: 


:29: 


:59. 


.5 


23{22 


5-75{I;gg 


3.o{|;| 


24{|I 


5-75{| : gg 


3-0{| 


4 





.026 




1 


SD 






1237-0151 






12: 


37: 


01. 


70 


-01: 


:51 


:23. 


.0 


23{| 


4.75{t;00 


3 .9{f:i 


25{|1 


4.75{f : oo 


3.9{| 


2 





.028 




1 






pWD 


1247-0039 


PG 1244- 


-004 


12: 


17: 


06. 


79 


-00: 


:39: 


:25. 


.8 


33{32 




1.8{f 3 s° 


39{f 7 


o-ooia 


2.i{«; 


.0 





.033 


WD 


3 


SD 






1300+0045 


PG 1257+010 


13: 


00: 


25. 


52 


+00: 


:45: 


:30 


.1 


29{i 


4.75{i;00 


i-4{?:i 


29{1 




i.2{i 


.9 





.026 




1 


SD 




pWD 


1300+0057 


HE 1258+0113 


13: 


00: 


59. 


21 


+00: 


:57. 


:11. 


.8 


3oqi 


3.50{|;il 


i-7{l:l 


32{|? 


3-50{|:I| 


1.8{1 


5 


0. 


.025 


WD 


1 


SD 


8 


SD 


1312+2245 






13: 


12: 


42. 


62 


+22: 


:45: 


:0d 


.2 


26{| 


3.00{° : °° 


3-7{I:I 


29{|g 


aoo{ ll° 


2.1{| s 


.1 


0. 


.013 




3 








1315+0245 






13: 


15: 


12. 


39 


+02: 


:45: 


31. 


.7 


33{32 


o nc 1-3.50 
J - zo l3.00 


o.9{J : g 


36{I 


3-25{ 3 oo 


0.9{i 







.023 




1 






pWD 


1316+0348 


PG 13144 


-041 


13: 


16: 


38. 


.48 


+03: 


:48: 


:18. 


.5 


29{lg 


4.75{t : oo 


l-5{fi 


29{i 


d.SOft?! 


1.3{1 


!) 





.030 




1 


SD 






1316+0739 






13: 


16: 


33 


59 


+07: 


:39: 


:41 


.3 


29{i 




6.o{^i 


29{|g 


6-50{| : | 


5.2{| 


4 


0. 


.026 




1 


SD 






1319-0141 






13: 


19: 


32. 


20 


-01: 


:41: 


:31. 


.2 


30{3i 


6-25{i:§8 


7.8{|| 


29{|| 


6.00{f;?l 


6.6{§ 


!) 


0. 


.024 




1 


SD 






1323+2615 






13: 


23: 


:57. 


28 


+26: 


:15: 


:02 


.5 


20{?^ 


5.00i1:|§ 


5 -8{3:6 


21{i 


4-25{| : gg 


3.3{| 


4 





.018 




1 






pWD 


1325+1212 


PG 13234 


-125 


13: 


25: 


:57. 


2d 


+ 12: 


:12: 


:21 


.3 


26{I 


5-75{|;gg 


2-i{?:l 


27{|I 


5-75{g : gg 


2.1{ 2 


3 
') 





.034 




1 








1326+0357 


PG 13234 


-042 


13: 


26: 


19. 


95 


+03: 


:57. 


:5d 


.4 


22{|? 


4.75{t;0° 


i-4{l:I 


25{| 


5-00{| : f5 


1.5{1 


() 


0. 


.025 


sdO 


1 


SD 






1328+3108 






13: 


28: 


:56. 


72 


+31: 


:08: 


:46. 







4-75{t:gg 


4.2{|:| 


38{ 37 


4.50{*:I| 


3.8{i 


" 


0. 


Oil 




d 


Galaxy 






1336+1126 


PG 13344 


-117 


13: 


36: 


53. 


99 


+ 11: 


:26: 


:05 


.4 


29{i 


4-25{|;gg 


Wtt 


30{3i 


4.25{| : | 


1.6{i 


.8 





.031 




1 








1341+0317 






13: 


41: 


22. 


97 


+03: 


:17: 


:51 


.6 


28{H 


6-50{|;75 


8.9{§ : g 


27{|1 


6-25{e:oo 


7.8{§ 


7 

') 





.024 


WD 


1 


SD 






1351+0234 






13: 


51: 


40. 


69 


+02: 


:34: 


:29. 


.2 


29{i 


4-25{l;§° 


2.5{f : | 


30{3i 


3-75{f|8 


2.4{| 


r> 


0. 


.027 


WD 


1 


SD 






1352+0910 






13: 


52: 


28. 


Id 


+09: 


:10: 


:39. 


.1 


29{i 


4-00{|;?I 


4-o{I:§ 


29{i 


4-oo{f : f| 


3.9{| 


3 

:) 


0. 


.028 


DA+M 


d 


cv 




WDMS 


1402+0725 


PG 13594 


-077 


Id: 


02: 


03. 


86 


+07: 


:25: 


39. 


.1 


26{i 


5.00{t;?| 


2-3{l:f 


2?{ 1 


5-OOftg 


2.3{| 


4 


0. 


.024 




1 


SD 






1402+3215 






Id: 


02: 


32. 


83 


+32: 


:15: 


:22 


.2 


22{|f 


6.25{g : §8 


i-9{?:l 


23(22 


6-25{e:oo 


1.9{f 


1 





.015 




1 


SD 






1421+0753 


KNBoo 
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